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ABSTRACT 

Data from the Fermz-LAT reveal two large gamma-ray bubbles, extending 50 degrees above and 
below the Galactic center, with a width of about 40 degrees in longitude. The gamma-ray emission 
associated with these bubbles has a significantly harder spectrum {dN/dE ~ E~'^) than the IC 
emission from electrons in the Galactic disk, or the gamma-rays produced by decay of pions from 
proton-ISM collisions. There is no significant spatial variation in the spectrum or gamma-ray intensity 
within the bubbles, or between the north and south bubbles. The bubbles are spatially correlated 
with the hard-spectrum microwave excess known as the WMAP haze; the edges of the bubbles also 
line up with features in the ROSATX-idiy maps at 1.5 — 2 keV. 

We argue that these Galactic gamma-ray bubbles were most likely created by some large episode of 
energy injection in the Galactic center, such as past accretion events onto the central massive black 
hole, or a nuclear starburst in the last ~ 10 Myr. Dark matter annihilation/decay seems unlikely 
to generate all the features of the bubbles and the associated signals in WMAP and ROSAT; the 
bubbles must be understood in order to use measurements of the diffuse gamma-ray emission in the 
inner Galaxy as a probe of dark matter physics. Study of the origin and evolution of the bubbles also 
has the potential to improve our understanding of recent energetic events in the inner Galaxy and the 
high-latitude cosmic ray population. 
Subject headings: gamma rays — ISM: jets and outflows — galaxies: active — galaxies: starburst 



1. INTRODUCTION 

The inner Milky Way is home to a massive black hole 
(MBH), surrounded by cl usters of youngj; stars and giant 
molecular clouds (see e.g. Morris fc Serabvnlll996[ for a 
review). The nuclear star cluster has a half-light radius 
of ~ 5pc. Although there are indications of past activity, 
the BH is quiescent today. 

Fe Ka echoes from molecular clouds around Sgr A* 
have been understo od as relics of activ i ty in the past 



1.1. Previous High-Energy Excesses 
Observations of gamma-ray emission in the inner 



Galaxy at £* < 1 GeV ^o back decade s to COS-B (IStroni 



19841: IStrong et al.ll987^ . and SAS-2 (IFichtel et al.lll975L 
Kniffen fc Fichtell Il98ll i (see iBloemenI Il989[ for a re- 
view). Later data from the EGRET experiment aboard 
the Compton Gamma-ray Obse rvatory extended to the 
high-energy side o f the tt^ bump (|Smialkowski et al.lll997 : 



Kovama et al. 



few h undred years (jSunvaev et al.l 119 ^ 

19961 ). On a longer timescale, one might expect relics of 
past activity in high energy CRs and hot gas, perhaps far 
off the disk. The most obvious observables would be e~ 
CR (visible in inverse Compton gammas and microwave 
synchrotron) and thermal emission (X-rays). 

This work presents a multiwavelength study of the in- 
ner Galaxy and identifies several large-scale (tens of de- 
grees) gamma-ray features, most notably 2 large (span- 
ning —50° < b < 50°) structures that we refer to as 
the ''Fermi bubbles". We suggest that these bubbles 
are associated with previously discovered structures in 
the X-rays and microwaves, and possibly with analogous 
smaller-scale structures visible in the FIR. 



Dixon et al.lll998f ). However, EGRET lacked the sensi- 
tivity and angular resolution to reveal the detailed struc- 
ture of gamma-ray emission toward the inner Galaxy. 
The Fermi Gamma-ray Space Telescope provides greatly 
improved data up to ~ 100 GeV, with sufficient angular 
resolution to map out interesting structures.^ 

At lower energies, the R OSAT All-Sky Survey at 1.5 
keV (jSnowden et al.l Il997l ) revealed a biconical X-ray 
structure over the inner tens of degrees around the Galac- 
tic center (GC), later interpreted as a superwind bub- 
ble (SWB) wi t h energetics of the order of IQ^^- ^^ ergs 
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( Sofud boOOat iBland-Hawthorn fc Cohenl l2003al ). On 
smaller scales, the Midcourse Space Experiment com- 
bined with IRAS data also confirms the existence of a 
limb-brightened bipolar structure, the so called Galactic 
center lobe (GCL), with orig in at the GC on the de- 
gree scale (see e.g. lLawll2010[ for a summary of multi- 
wavelength observations of GCL). The inferred energy 
injection of both these bipolar structures, despite their 

^ See http://fermi.gsfc.nasa.gov/ssc/data/ 



different scales, is ~ 10^^ ^^ ergs, with an estimated age 
of - 10^ yr for the GCL and - 10^ yr for the SWB. 
Several Galactic Center Shells (GCS), tens of pc in size, 
have been fo und with total energy of order ~ 10^^ ergs 
(|Sofuell2QQ3l ). These shells and filaments are claimed to 
originate from one or more episodes of rapid energy re- 
lease. 

1.2. Microwave Excess: the WMAP Haze 

Beyond direct evidence of shell structures, microwave 
observations also provide intriguing indications of energy 
release toward the GC. 

At tens of GHz, the Wilkinson Microwave Anisotropy 
Probe ( WMAP) ^ provides sensitive degree resolution full 
sky maps of diffuse microwave emission. By subtracting 
templates including Galactic Ha, Haslam 408 MHz soft 
synchrotron, and thermal dust emission to remove the 
different known emission mechanisms in these maps, a 
microwave residual excess (named "the microwave haze" ) 
with spherical (non-disklike) morphology about ^ 4 kpc 
in radius toward the GC (visible up to a t least \b\ ^ 30°) 
has been recognized (jFinkbeinerl l2004al ) . It has a spec- 
trum of about Iiy ~ i/~^-^, harder than typical syn- 
chrotron, but softer than free- free. The microwave haze 
was later interpreted as synchrotron emission from a hard 
spectrum of e~ cosmic rays. Other hypotheses such as 
free- free, spinning dust, or thermal dust ha ve failed to 
explain its morphology, spectruni , or both (jFinkbeiner 
2QQ4b[ bobler fc Finkbeinerl [2OO8I ) . However, the most 
recent WMAP 7- year data have not detected the haze 
polarization pred icted by the synchrotron hypothesis 
(jGold et al.ll2010[ ). implying either heavily tangled mag- 
netic fields, or an alternative emission mechanism. With 
that caveat in mind, we will assume the WMAP haze is 
synchrotron and consider the implications. 

1.3. A Hard Electron CR Spectrum 

A simple model, in which the electron CRs that form 
the haze have diffused from supernova shocks in the disk, 
cannot fully explain the data for standard diffusion as- 
sumptions. The 23 — 33 GHz spectrum of the haze 
synchrotron is as hard as that generated from shocks, 
and it seems extremely unlikely that these electrons 
can diffuse several kpc from the disk without signifi- 
cant softening of the spectrum. The synchrotron cooling 
timescale for cosmic ray electrons emitting at frequency 
V is Tsvn ^ 10^ ^,r^^ ^ CH- y^' where 5ioo = 5/100/iG 
(jThompson et al.ll2006[ ). Besides the hard spectrum, it is 
difficult to form the distinctly non-disklike morphology 
of the haze with any population of sources concentrated 
in the disk (as is believed to be true of supernovae). 

The presence of a distinct component of diffuse hard 
e~ CR far off the plane is intriguing in itself, and has 
motivated proposals where the haze is generated by 

^ http://inap.gsfc.nasa.gov/ 



pulsars, other astroph ysical processes, or the annihila- 



i Kaplinghat et a 



iMalvshev et al.l 



tion o f dark matte r (IHooper et al.l 120071: ICholis et al. 
^2009a! ; Izhang et al.l [2OO9I: iHarding fc AbazaiianI I2OIO : 



1 I2OO9I : iMcQuinn fc Zaldarriagal I2OIO : 



201 Of ) . Other indications of excess elec- 



tronic activity in the Milky Way may be found in re- 
cent measurements of local electron and positron cos- 
mic rays. The ATIC, Fermi and H.E.S.S experiments 
have observed a h ardening in the e^ + e~ spectrum 
at 20 - 1000 GeV dchang et al.l [20081 lAharonian et al' 



2009l : lAbdo et al.ll2009[ ). with an apparent steepening at 
~ 1 TeV, and the PAMELA experiment has measured a 
rising positron fraction above 10 GeV. Taken together, 
these measurements imply a new source of hard electrons 
and positrons, which may be related to the WMAP haze. 
The coexistence of ROSAT X-ray bipolar features and 
the WMAP haze toward the inner Galaxy also suggests 
the interesting possibility of a common physical origin 
for these signals. 

1.4. Inverse Compton Excess from Fermi-LAT 

Fortunately, if the WMAP haze is synchrotron radi- 
ation from a hard electron population located around 
the GC, the same CRs would also produce IC scattered 
gammas, allowing an independent probe of the CR pop- 
ulation. IC photons provide valuable complementary in- 
formation about the spatial distribution of the e~ CR 
(given a model for the ISRF), which in turn can con- 
strain hypotheses about their origin. The unprecedented 
angular resolution and sensitivity of Fermi-LAT allows 
us to probe the gamma-ray counterpart to the microwave 
haze in detail for the first time. 

Previous work employing the first year Fermi-LAT 
data isolated a spectrally hard "gamma-ray haze" with 
similar morphology to the WMAP microwave haze 
(JDobler et al.l l2010[ ). In this work, we show that the 
Fermz-LAT sky maps constructed from 1.6 yr data (600 
days) reveal two large gamma-ray lobes, extending 50 
degrees above and below the GC, with a width of about 
40 degrees in longitude. These two "bubble" -like struc- 
tures have relatively sharp edges and are symmetric with 
respect to the galactic plane and the minor axis of the 
galactic disk. The gamma-ray signal reveals similar mor- 
phology to the WMAP haze, and is also suggestive of a 
common origin with features in the i^O^^dT X-ray maps 
at 1.5 keV towards the GC. 

As we will discuss, the sharp edges, bilobular shape, 
and apparent centering on the GC of these structures 
suggest that they were created by some large episode of 
energy injection in the GC, such as a past accretion onto 
the central black hole, or a nuclear starburst in the last 
~10 Myr. It is well known that the GC hosts a massive 
black hole and massiv e clusters of recently formed stars 
( Paumard et al.ll2006[ ). Either of these could potentially 
provide the necessary energy injection by driving large- 
scale galactic winds or producing energetic jets; we will 



outline some of the advantages and disadvantages of each 
scenario. 

1.5. Structure of This Paper 

In Section [2] we briefly review the Fermi-LAT data and 
our data analysis procedure. In Section [3l we show the 
1.6 yr Fermi data maps, reveal the Fermi bubble features 
and show that they are robust when different models for 
the expected Galactic diffuse emission are subtracted. 
We characterize the morphology of the bubbles in some 
detail and employ regression template fitting to reveal 
a hard, spatially uniform spectrum for the gamma-ray 
emission associated with the bubbles. In Section [H we 
show that features in the ROSAT soft X-rays and the 
WMAP microwave haze are spatially correlated with the 
Fermi bubbles, and the WMAP haze and Fermi bub- 
bles are consistent with being produced from the same 
electron CR population (by synchrotron and IC respec- 
tively). Section [5] presents our understanding based on 
the analysis in Section [3] and Section [H Section [6] dis- 
cusses possible scenarios to produce the Fermi gamma- 
ray bubbles. Section [71 focuses on the origin of the elec- 
tron CRs and the challenges in explaining the spectral 
and spatial profiles of the gamma-ray emission from the 
bubbles. We discuss the implications of the Fermi bub- 
bles for several topics of interest in Section [51 and give 
our conclusions and suggest future work in Section [9l 



2. FERMI DATA AND MAP MAKING 
The Large Area Telescope (^LAT; 



, see 
Gehrels & MichelsonI Il999l : lAtwood et aJ] I2QQ9I : as 
well as the Fermi homepage^) is a pair-conversion 
telescope consisting of 16 modules of tungsten and 
silicon-strip trackers, on top of a calorimeter with a 
thickness of 7 radiation lengths. It has a scintillating 
anti-coincidence detector that covers the tracker array, 
and a programmable trigger and data acquisition sys- 
tem. The LAT provides a wide field of view, and covers 
the energy range from about 30 MeV to 300 GeV. 

The spacecraft occupies a low Earth orbit with an in- 
clination of 25.6°. The field of view is so wide that 
the entire sky may be covered in two orbits by rocking 
the spacecraft north of zenith on one orbit and south of 
zenith on the next. This scan strategy exposes the LAT 
to atmospheric gammas at high zenith angles. We make 
use of only events designated "Class 3" (P6_V3 diffuse 
class) by the LAT pipeline with a zenith angle cut of 
105°. 

The events are binned into a full sky map using 
HEALPix, a convenient iso-latitude equal-area full-sky 
pixelization widely used in the CMB community.^ Spher- 

^ http://fermi.gsfc.nasa.gov/ 

^ HEALPix software and documentation can be found at 
http://healpix.jpl.nasa.gov, and the IDL routines used in 
this analysis are available as part of the IDLUTILS product at 
http : //sdssSdata . Ibl . gov/sof tware/idlutils. 



ical harmonic smoothing is straightforward in this pix- 
elization, and we smooth each map by the appropriate 
kernel to obtain a Gaussian PSF of 2° FWHM. Because 
the PSF of the initial map must be smaller than this, 
at energies below 1 GeV we use only front-converting 
events (which have a smaller PSF). A larger smooth- 
ing scale would help improve S/N, but a relatively small 
smoothing scale is necessary to see sharp features (such 
as the bubble edges). Furthermore, for the comparisons 
and linear combination analysis described in e.g. §3.1.31 
it is necessary to smooth the maps at each energy to 
a common PSF. We generate maps from 1 — 4°, and 
find that a FWHM of 2 ° works well for our purposes. 
See IPobler et al.l (|2Q1Q[ ) for details on map construc- 
tion, smoothing, masking, and for instructions on how 
to download the maps. 

Our current gamma-ray maps (v2_3) constructed from 
the Fermi data have greater signal/noise compared 
to the previously released vl_0 maps. They contain 
photon events from mission times 239557417.494176 to 
291965661.204593, for about 606 days or 1.66 years of 
data (rather than 1 year). We refer to th ese as the "1.6 



year maps." As in iDobler et al.l ()2010l ). we construct 



maps of front-converting and back-converting events sep- 
arately, smooth to a common PSF, and then combine 
them. The point source subtraction has improved: in- 
stead of interpolating over every source in the 3-month 
catalog, we use the 1-year catalog,^ and subtract each 
point source from the maps in each energy bin, using es- 
timates of the PSF from the Fermi science tools. ^ For 
the 200 brightest and 200 most variable sources, the sub- 
traction is noticeably imperfect, so we interpolate over 
the core of the PSF after subtracting the best estimate. 
We take care to expand the mask for very bright sources 
(Geminga, 3C 454.3, and LAT PSR J1836+5925). The 
resulting map is appropriate for diffuse work at \h\ > 3°. 
At \b\ < 3° the maps are severely compromised by the 
poor subtraction and interpolation over a large number 
of point sources. Further details of the map process - 



ing may be found in Appendix B of lDobler et al 



(|2010h . 

The v2_3 maps used in this work and color versions of 
the map are available for download. ^^ 

3. FERMI BUBBLES 

3.1. Diffuse Galactic Emission Models 

At low (~1 GeV) energies, and close to the Galactic 
plane, the gamma-rays observed by Fermi are dominated 
by photons from the decay of tt^ particles, produced by 
the collisions of CR protons with ambient ionized gas 
and dust in the ISM. Collisions of high energy CR elec- 
trons with the ISM (primarily protons, but also heavier 
nuclei) produce bremsstrahlung radiation. The Fermi 

^ Available from http://fermi.gsfc.nasa.gov/ssc/data 

^ See http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/ 

^^ Available at http : //f ermi . skymaps . info 
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Fig. 1. — All-sky Fermi-hAT 1.6 year maps in 4 energy bins. Point sources have been subtracted, and large sources, including the inner 
disk (-2° < 6 < 2°, -60° <i < 60°), have been masked. 



all-sky gamma-ray maps in different energy bands are 
shown in Figure [H In order to uncover the Fermi bubble 
features better, significant tt^ emission, bremsstrahlung, 
and IC emission from the Galactic disk must be removed. 
We take three approaches for the foreground removal. 
One is to use the Fermi Diffuse Galactic Model provided 
by the Fermi team^^( §3.1.ip . The second approach em- 
ploys a linear combination of templates of known emis- 
sion mechanisms ( ^3.1.2p . using existing maps from mul- 
tiwavelength observations and/or constructed geometric 
templates. The third approach is taking advantage of 
the lower energy band 0.5 — 1.0 GeV Fermi map to form 
a template of a diffusion emission model ( §3.1.3p . 

3.1.1. Fermi Diffuse Galactic Model 

The Fermi diffuse Galactic model^^ is a comprehen- 
sive model of Galactic gamma-ray emission from the 
ISM, and serves as a background estimate for point 
source removal. This model is based on template fits 
to the gamma-ray data, and includes an IC compo- 
nent generated by the GALPROP cosmic ray propagation 
code. GALPROP calculates the steady state solution to 
the diffusion-energy-loss equation, given the 3D gas dis- 
tribution, interstellar radiation field, B-field model, CR 
diffusion assumptions, and many other input par a meter s 
( Strong fc Moskalenkol Il999l : IStrong et all 120091 120071 ). 



The model is constrained by gamma-ray and microwave 
observations, locally measured CR spectra, etc. By us- 
ing a well motivated physical model, one can solve for 
the spectral and spatial dependence of the injection func- 
tion, i.e. the e~ and p CR primary source spectra, as a 
function of position and energy. The diffuse model is 
the key connection between the input assumptions and 
the observables, and is essential for interpretation of the 
Fermi-LAT data. It is important to make it as complete 
as possible. 

In this model, the tt^ emission is modeled with maps 
of interstellar gas: H i fro m the Leiden/ Argent ine /Bonn 
(LAB) Galactic Survey (|Kalberla etld] l2005h a n^CO 
from the CfA composite CO survey (JDame et al.ll200l[ ). 
Because the tt^ emission is a function of both the gas den- 
sity and the proton CR density, which varies with Galac- 
tocentric radius, it is desirable to allow the emissivity 
of the gas to vary. Both the H i and CO surveys con- 
tain velocity information, which allows separation into 
six Galactocentric annuli (rings) with boundaries at 4.0, 
5.5, 7.0, 10.0, 16.5, and 50 kpc. The spectrum of each is 
allowed to float, with the constraint that the sum of the 
rings along each line of sight approximates the observed 
signal. This freedom also allows for varying amounts 
of bremsstrahlung (with varying spectrum) which also 
scales with the ISM density. The contribution from IC is 
modeled with GALPROP as described above, and included 



^^ See http : //f ermi . gsf c . nasa. gov/ssc/data/access/lat/Backgroun!(iM6'<M.^.l&finP-^ 
^"^ Available from http://fermi.gsfc.nasa.gov/ssc/data. The 
version of the diffuse model we use is gll.iem.v02. ^ description 
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Fig. 2. — All-sky residual maps after subtracting the Fermi diffuse Galactic model from the LAT 1.6 year maps in 4 energy bins (see 
^3.1.1|) . Two bubble structures extending to 6 ± 50° appear above and below the GC, symmetric about the Galactic plane. 



This procedure provides a diffuse model that faithfuhy 
reproduces most of the features of the diffuse Galactic 
emission. One shortc oming is the existence of "dark gas" 
( Grenier et alJl2QQ5l ). clouds with gamma-ray emission 
that do not appear in the Hi and CO surveys. T hese 
features are seen in dust maps (jSchlegel et al.lll998[ ) and 
may simply be molecular H clouds under abundant in CO. 

The Fermi diffuse model is primarily intended as a 
background for point source detection, and comes with a 
number of caveats. However these caveats apply mainly 
near the Galactic plane, and at £^ > 50 GeV. It is nev- 
ertheless useful for qualitatively revealing features in the 
diffuse emission at high latitude. In Figure [21 we show 
the residual maps after subtracting the Fermi diffuse 
Galactic model in different energy bins. A double- lobed 
bubble structure is clearly revealed, with similar mor- 
phology in the different energy bins. We note that the 
bubble is neither limb brightened nor centrally bright- 
ened, consistent with a flat projected intensity distribu- 
tion. 

3.1.2. Simple Template- Based Diffuse Galactic Model 

Since the dominant foreground gamma-rays originate 
from TT^ gammas produced by CR protons interact- 
ing with the ISM, the resulting gamma-ray distribution 
should be morphologically correlated with other maps 
of spatial tracers of the ISM. A good candidate is the 
Schlegel, Finkbeiner, & Davis (SFD) map of Galactic 

ferini.gsfc.nasa.gov/ssc/data/access/lat/ring_for_FSSC_final4.pdf 



dust, based on 100 jjlvh far IR data (jSchlegel et al.lll998[ ). 
The TT^/bremsstrahlung gamma-ray intensity is propor- 
tional to the ISM density x the CR proton/electron den- 
sity integrated along the line of sight. As long as the 
CR proton/electron spectrum and density are approxi- 
mately spatially uniform, the ISM column density is a 
good tracer of 7r^/bremsstrahlung emission. The dust 
map has some advantages over gas maps: there are no 
probl ems with self absorp tion, no concerns about "dark 
gas" (jGrenier et al.ll2005[ ). and the SFD dust map has 
sufficient spatial resolution (SFD has spatial resolution 
of 6', and LAB is 36'). On the other hand, SFD con- 
tains no velocity information, so it is impossible to break 
the map into Galactocentric rings. Nevertheless, it is in- 
structive to employ the SFD map to build a very simple 
foreground model. The goal is to remove foregrounds in a 
fashion that reveals the underlying structure with as few 
physical assumptions as possible. We will compare the 
resulting residuals using this simple diffuse model with 
those using the Fermi diffuse Galactic model. 

As an example, we reveal the Fermi bubble structure 
from 1 — 5 GeV Fermi-LAT 1.6 yr data in Figure [3l We 
use the SFD dust map as a template of the tt^ gamma 
foreground. The correlation between Fermi and SFD 
dust is striking, and the most obvious features are re- 
moved by this subtraction {top row in Figure [3]). This 
step makes the bubbles above and below the GC easily 
visible. The revealed bubbles are not aligned with any 
structures in the dust map, and cannot plausibly be an 
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Fig. 3. — Template decomposition of the Fermi-L AT 1.6 year 1 — 5 GeV map (see ^3.1.2|) . Top left: Point source subtracted 1 — 5 
GeV map, and large sources, including the inner disk (—2° < 6 < 2°,— 60° < i < 60°), have been masked. Top middle: The 1 — 5 
GeV map minus SFD dust map {top right panel) which is used as a template of tt*^ gammas. Middle row: The left panel is the same as 
the top middle panel but stretched 2x harder. The m^iddle panel subtracts a simple geometric disk template (shown in the right panel), 
representing mostly inverse Compton emission, to reveal features close to the Galactic center. Two large bubbles are apparent (spanning 
—50° < 6 < 50°). Bottom, row: The left panel is the same as the middle panel of the second row. Finally we subtract a simple bubble 
template {right panel), with a shape derived from the edges visible in the maps, and uniform projected intensity. After subtracting the 
bubble template, the two bubbles features have nearly vanished {bottom, m^iddle panel), indicating a nearly flat intensity for the Ferm^i 
bubbles. 
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Fig. 4. — Full sky residual maps after subtracting the SFD dust and disk templates from the Fermi-LAT 1.6 year gamma-ray maps in 
two energy bins. Point sources are subtracted, and large sources, including the inner disk (—2° < b < 2°,— 60° < i < 60°), have been 
masked. Two large bubbles are seen (spanning —50° < b < 50°) in both cases. Right panels: Apparent Fermi bubble features marked 
in color lines, overplotted on the maps displayed in the left panels. Green dashed circles above and below the Galactic plane indicate the 
approximate edges of the north and south Fermi bubbles respectively. Two blue dashed arcs mark the inner (dimmer) and outer (brighter) 
edges of the northern arc - a feature in the northern sky outside the north bubble. The red dotted line approximately marks the edge of 
Loop I. The purple dot-dashed line indicates a tentatively identified "donut" structure. 



artifact of that subtraction. 

Next, a simple disk model is subtracted (Figure [3l mid- 
dle row). The purpose of this subtraction is to reveal the 
structure deeper into the plane, and allow a harder color 
stretch. The functional form is (csc|6|) — 1 in latitude 
and a Gaussian {(j£ = 30°) in longitude. The disk model 
mostly removes the IC gamma-rays produced by cosmic 
ray electrons interacting with the ISRF including CMB, 
infrared, and optical photons; as discussed previously, 
such electrons are thought to be mostly injected in the 
Galactic disk by supernova shock acceleration before dif- 
fusing outward. 

Finally, we fit a simple double-lobed geometric bub- 
ble model with flat gamma-ray intensity to the data, to 
remove the remaining large-scale residuals towards the 
GC (Figure [31 bottom row). In this model, we identify 
the approximate edges of the two bubble-like structures 
towards the GC in the bottom left panel (shown with 
dashed green line in right panels of Figure [4]) . We then 
fill the identified double-lobed bubble structure with uni- 
form gamma ray intensity, as a template for the ''^ Fermi 
bubbles" {bottom right panel of Figure [3]). If the Fermi 
bubbles constitute the projection of a three dimensional 
two-bubble structure symmetric to the Galactic plane 
and the minor axis of the Galactic disk, taking the dis- 
tance to the GC Rq = 8.5 kpc, the bubble centers are 



approximately 10 kpc away from us and 5 kpc above and 
below the Galactic center, extending up to roughly 10 
kpc as the most distant edge from GC has \b\ ~ 50°. 
No structures like this appear in GALPROP models, and in 
fact GALPROP is often run with a box-height smaller than 
this. Because the structures are so well centered on the 
GC, they are unlikely to be local. 

In Figure [H we show the full sky residual maps at 1 — 5 
GeV and 5 — 50 GeV after subtracting the SFD dust and 
the disk model to best reveal the Fermi bubble features. 
Although photon Poisson noise is much greater in the 
5 — 50 GeV map, we identify a Fermi bubble structure 
morphologically similar to the structure in the 1 — 5 GeV 
map, present both above and below the Galactic plane. 

In Figure[5l we show the full sky maps at 1—5 GeV with 
the zenithal equal area (ZEA) projection with respect to 
both north pole and south pole. We found no interesting 
features appear near the poles. 

3.1.3. Low Energy Fermi Map as a Diffuse Galactic Model 

In Figure El we show the 0.5 - 1 GeV and 2-50 GeV 
residual maps after subtracting only the SFD dust map 
as a template of foreground tt^ gammas. The residual 
maps should be dominated by IC emission from CR elec- 
trons interacting with the ISRF. We use the 0.5 — 1 GeV 
maps as a template of IC emission from high energy elec- 
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Fig. 5. — The zenithal equal area (ZEA) projection with respect to 
the 1 — 5 GeV energy band before {left panels) and after (right panels) 
Fermi bubbles. 

trons scattering starlight^ and subtract the template from 
higher energy maps (the lower panels of Figure [6|). The 
Fermi bubble structures are clearly revealed. We thus 
conclude that the Fermi bubbles are mostly from high 
energy electron CRs IC scattering on CMB photons, and 
IR photons at higher energies (see ^ for more discus- 
sion). By comparing the Fermi diffuse Galactic model 
subtraction (Figure [2]) and our simple template model 
subtraction (Figure [4]), we find that the bubble structures 
are robust to quite different foreground subtractions. It 
is difficult to see how such emission could arise - espe- 
cially with sharp edges Section 13.2.11 - as an artifact of 
these subtractions. 

3.2. Fermi Bubbles: Morphology 
3.2.1. Morphological Features 

In the right panels of FigureHl we illustrate the edges of 
the Fermi bubbles and some other features. We find that 
the Fermi bubbles have distinct sh arp ed^es, rath er than 
smoothly falling off as modeled in iDobler et al.l (|2QlQi ). 



both north pole (upper panels) and south pole (lower panels) for 
subtracting the SFD dust and simple disk templates to reveal the 

Besides the two bubbles symmetric with respect to the 
Galactic plane, we find one giant northern arc that em- 
braces half of the north bubble, that extends from the 
disk up to 6 ~ 50, with i ranging from roughly —40° 
to 0°. It has a brighter and sharper outer edge in the 
1 — 5 GeV map. On a even larger scale, we identify a 
fainter structure extended up to 6 '-«^ 80°, with £ ranging 
from roughly —80° to 50°. We will show in §4lthat this 
large extended structure corresponds to the North Polar 
Spur emission associated with Loo p I (as seen for ex - 
ample in the Haslam 408 MHz map lHaslam et al.lll982[ ). 
We will discuss the possible relation of the Fermi bubble, 
the northern arc^ and the Loop I feature in ^ In the 
1 — 5 GeV map, we also identify a "donut-like" struc- 
ture in the south sky with b ranging from roughly —35° 
to 0° and £ from roughly 0° to 40°. The coordinates of 
the Fermi bubble edges, northern arc^ Loop /, and the 
"donut" feature identified from the 1 — 5 GeV map are 
listed in Table [B 
In Figure [71 we compare the Fermi bubble morphology 
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Fig. 6. — Top left: Full sky Fermi-L AT 1.6 year 0.5 — 1.0 GeV map subtracts the SFD dust map as a template of 7r° gammas. Top 
right: The same as top left panel, but for energy range 2 — 50 GeV (note the different gray scale for the two panels). Bottom left: The 
2 — 10 GeV Fermi gamma-ray map subtracting the top left 0.5 — 1.0 GeV residual map which is used as a template of ICS of starlight. 
Bottom, right: The same as bottom, left panel but for 10 — 50 GeV map subtracting the top left 0.5 — 1.0 GeV residual map. The Fermi 
bubble structures are better revealed after subtracting the lower energy 0.5 — 1.0 GeV residual map with extended disk-like emission. 

in different energy bins. We show the difference of the 
1 — 2 and 2 — 5 GeV residual maps in the upper panels] 
each residual map is the result of subtracting the SFD 
dust map and the simple disk model to best reveal the 
Fermi bubbles. The difference maps between the 1 — 5 
and 5 — 50 GeV maps are shown in the lower panels. The 
bubble features almost disappear in the difference maps, 
indicating that different parts of the Fermi bubbles have 
similar spectra. 

To study the sharp edges of the bubbles at high lati- 
tude more carefully, we examine the (projected) intensity 
profiles along arcs of great circles passing through the 
estimated centers of the north and south bubbles, and 
intersecting the bubble edge (as defined in Figure [4]) at 
|6| > 28°. Along each such ray, we define the intersection 
of the arc with the bubble edge to be the origin of the 
coordinate system; we then perform an inverse- variance- 
weighted average of the intensity profile along the rays 
(as a function of distance from the bubble edge). We 
subtract a constant offset from the profile along each 
ray, prior to averaging the rays together, to minimize 
aliasing of point sources onto the averaged profile, and 
then add the averaged offset back in at the end. The 
inverse variance for each data point is obtained from the 
Poisson errors in the original photon data, prior to any 
subtraction of point sources or templates (however, the 
smoothing of the map is taken into account). When the 



rays are averaged together, the naive inverse variance in 
the result is multiplied by a factor of the annulus ra- 
dius (for the points being averaged together) divided by 
47rcr^, where cr is the la value of the PSF, and the an- 
nulus width is taken to be 1° (the spacing between the 
points along the rays; this is comparable to the smooth- 
ing scale, so there may still be unaccounted-for correla- 
tions between the displayed errors); this is done to take 
into account that the number of independent measure- 
ments being sampled by the rays can be far less than 
the number of rays, especially close to the center of the 
bubbles. This procedure is repeated for all the stages of 
the template subtraction, using the simple disk template 
for inverse Compton scattering (ICS) for illustration (our 
conclusions do not depend on this choice). 

The results are shown in Figure [8] for the averaged 
(1 — 2) + (2 — 5) GeV maps, and the averaged (5 — 10) + 
(10 — 20) GeV maps. In both energy ranges the edges 
are clearly visible; in the south, this is true even before 
any templates are subtracted. The intensity profile of 
the north bubble is strikingly similar to profile of the 
south bubble. For both of the north and south bubbles, 
no significant edge-brightening or limb-brightening of the 
bubbles is apparent from the profiles, the flux is fairly 
uniform inside the bubbles. 

In Figure [H we plot the intensity profile as a function 
of latitude from the south to the north pole. We con- 
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Fig. 7. — Comparison of the Fermi bubbles between different energy bins. Top left: Subtraction of the 2 — 5 GeV residual map from 
the 1 — 2 GeV residual map; each residual map is constructed from the data by regressing out the SFD dust and disk templates to best 
reveal the Fermi bubbles. The difference map is consistent with Poisson noise away from the masked region, and the bubble features can 
hardly be recognized, indicating that different spatial regions of the Ferm,i bubbles have the same spectrum. Top right: The same map as 
left panel, but with the Ferm.i bubble features overplotted for comparison. The marked features are the same as those plotted in Figure 
m and are listed in Table [T] Bottom, row: Same as the upper panels, but subtracting the 5 — 50 GeV residual map from the 1 — 5 GeV 
residual map. 
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Fig. 8. — Intensity as a function of radial distance from the bubble edge, averaged over great circle arcs intersecting the bubble center 
and lying at |6| > 28°. Results are shown for (left) the southern bubble, and (right) the northern bubble, for (top) the averaged 1 — 2 
and 2 — 5 GeV maps, and (bottom) the averaged 5 — 10 and 10 — 20 GeV maps. Different lines show the results at different stages of the 
template regression procedure and the corresponding errors are plotted (see text for an outline of the error analysis). 



struct great circle arcs perpendicular to the / = great 
circle, extending 10° in each direction (east and west), 
and average the emission over each such arc. The flatness 
of the bubbles with latitude (except possibly close to the 
Galactic plane), and the sharp edges at high latitude, are 
also apparent here. 

We note that the flat intensity of the bubbles is strik- 
ing. As we show in Figure [TOl if we assume that the 
Fermi bubbles are projected from spherically symmetric 



three-dimensional bubbles centered above and below the 
GC, a non-trivial emissivity distribution in the bubble 
interior is required to produce a flat projected intensity 
distribution {upper left panel of Figure [TQ|) . If the "bub- 
bles" originate from IC scattering, this suggests a rather 
non-uniform density distribution for the electron CRs, 
which - combined with a nearly uniform spectral index 
- presents challenges for many models for the electron 
injection. The expected intensity profile for a shock gen- 
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Fig. 9. — Intensity averaged over the central 20 degrees in longitude, as a function of latitude, for (left) the averaged 1 — 2 and 2 — 5 GeV 
maps, and (right) the averaged 5-10 and 10-20 GeV maps. We construct great circle arcs perpendicular to the I = great circle, extending 
10° in each direction (east and west), and average the emission over each such arc; the "6" label corresponding to each arc, and the x-axis 
of the plot, refers to the value of 6 at / = 0. Different lines show the results at different stages of the template subtraction process. The 
large oversubtraction at 6 ~ 15° in the north, especially pronounced in the low-energy data, is associated with a bright feature in the SFD 
dust map. 

erated bubble with a compressed gas shell is shown in 
the upper right panel; it is noticeably limb-brightened, 
in contrast to observations. 



3.2.2. Spectrum of Gamma-ray Emission 

We now attempt to estimate the spectrum of the 
gamma rays associated with the Fermi bubbles, and the 
spatial variation of the spectrum. In order to reveal 
the hardness of the spectrum of the Fermi bubbles, and 
quantitatively study the intensity flatness of the bub- 
ble interiors, we do a careful regression template fitting. 
First, we maximize the Poisson likelihood of a simple dif- 
fuse emission model involving 5-templates (see §3.1.2p . In 
this model, we include the SFD dust map (Figure [31 right 
panel of the top row) as a tracer of tt^ emission which is 
dominant (or nearly so) at most energies on the disk and 
significant even at high latitudes, the simple disk model 
(Figure [3l right panel of the second row)^ the bubble 
template (Figure [3l right panel of the bottom row), the 
Loop I template (see e.g. Figure [TT| left panel of the top 
row)^ and a uniform background as templates to weight 
the Fermi data properly. 

The significant isotropic background is due to extra- 
galactic emission and charged particle contamination, in- 
cluding heavy nuclei at high energies. The Fermi collab- 
oration has measured the extragalactic diffuse emission 
using additional cuts to r educe charged particle contami- 
. Abdo et all (|2QlQh : below - 20 GeV the isotropic 



nation 



than the extragalactic diffuse emission, but has a similar 
spectral slope. At energies above ~ 20 GeV the isotropic 
contribution becomes much harder, which we attribute 
to charged particle contamination. 

For each set of model parameters, we compute the Pois- 
son log likehhood. 



In/: 






hlniii - iii - ln(/ci!). 



(1) 



contribution in our fits is roughly a factor of 2 larger 



where jii is the synthetic map (i.e., linear combination of 
templates) at pixel i, and k is the map of observed data. 
The last term is a function of only the observed maps. 
The 1(7 Gaussian error is calculated from the likelihood 
by Aln>C = 1/2. The error bars are simply the square 
root of the diago nals of the coyarianc e matrix. We refer 
to Appendix B of lPobler et al.l (|2010l ) for more details of 
the likelihood analysis. Maps of the models constructed 
from linear combinations of these five templates, and the 
residual maps between the Fermi data and the combined 
templates at different energy bins, are shown in Figure 
[TTl In this fit, we mask out all pixels with Galactic lat- 
itude \h\ < 30° (the dashed black line in the residual 
maps). 

Template-correlated spectra for the 5-template fit are 
shown in Figure [121 The fitting is done with regions 
of |6| > 30°. For a template that has units (e.g., the 
SFD dust map is in Eb-v magnitudes) the correlation 
spectrum has obscure units (e.g. gamma-ray emission 
per magnitude). In such a case we multiply the correla- 
tion spectrum by the average SFD value in the bubble 
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Fig. 10. — Projected emission from four emissivity distributions to illustrate the qualitative features of each. Top left: A toy model of 
the Fermi bubbles with flat projected intensity. In this model, the volume emissivity (assuming the ISM is optically thin to gamma rays) 
is proportional to (R^ — r^)~^/^ for r < R, and zero otherwise, where R=3.5 kpc is the approximate radius of the two bubbles, and r is 
the distance to the center of the north or south bubble. Top right: A bubble model with compressed gas shells with a thickness of 0.5 kpc; 
the electron CR density in the shell is a factor of 5 higher than in the interior of the bubbles. For this model, a limb-brightened edge of 
the bubbles is clearly visible, a feature which is not seen in the Fermi data. Bottom, left: An illustrative toy model for the WMAP haze. 
The haze synchrotron emissivity depends on the electron CR density and the magnetic fleld; here we take B = Boe~^^^^ , where zq = 2 
kpc, and show the line of sight integral of this B-fleld through bubble volume. Even though the synchrotron emissivity is not simply the 
product of CR density times fleld strength, this panel suggests that the decreasing intensity of the WMAP haze at high latitudes is due 
to the decay of the Galactic magnetic fleld away from the Galactic plane. Bottom, right: A toy model for the ROSAT ^-ray features. The 
observed soft X-rays are limb-brightened and we assume all gas is uniformly distributed within a compressed shell, with no contribution 
from the interior, and X-ray emission is proportional to the gas density squared. The thickness of the shell is 1 kpc. 



region, defined by the bottom right panel of Figure [31 
masking out the \b\ < 30° region. For the uniform, Loop 
/, and bubble templates (including inner, outer, north, 
and south), no renormalization is done. These templates 
are simply ones and zeros (smoothed to the appropriate 
PSF), so the outer bubble spectrum is simply the spec- 
trum of the bubble shell template shown in Figure [151 
not the mean of this template over the whole bubble re- 
gion. The normalization factors for different templates 
are listed in Table [21 

In Figure [121 we show spectra for tt^ emission, 
bremsstrahlung and inverse Compton scattering calcu- 
lated using a sample GALPROP model (tuned to match 
locally measured protons and anti-protons as well as lo- 
cally measured electrons at ~ 20 — 30 GeV), as an indi- 
cation of the expected spectral shapes. The spectra for 



the SFD and the simple disk template reasonably match 
the model expectations. The dust map mostly traces 
the TT^ emission, and the simple disk model resembles 
a combination of IC and bremsstrahlung emission. The 
spectrum for emission correlated with the Fermi bubbles 
is clearly significantly harder than either of these com- 
ponents, consistent with a flat spectrum in E'^dN/dE. 
This fact coupled with the distinct spatial morphology 
of the Fermi bubbles indicates that the IC bubbles are 
generated by a separate electron component. We also 
note that the spectrum of the bubble template falls off 
significantly at energy <1 GeV. This feature is robust 
with respect to the choice of templates. The fitting co- 
efficients and corresponding errors of each template are 
listed in Table [3l We will discuss some implications of 
the falling spectrum in ^4.31 
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To demonstrate the robustness of the spectrum we have 
derived for the Fermi bubbles, we make use of the Fermi 
0.5 — 1 GeV residual map (after subtracting the SFD dust 
map to largely remove the tt^ gammas) as a template of 
IC emission, and perform a 4-template fit ( §3.1.3p . These 
gamma rays mostly originate from IC scattering of a rel- 
atively soft population of electrons in the disk, but might 
also contain gammas from IC scattering on starlight by 
a latitudinally extended electron population. We use the 
SFD dust map as a template for tt^ gammas as previ- 
ously, and include the uniform background and the bub- 
ble template as in the previous 5-template fit. The fitting 
is done with regions of |6| > 30°. For the SFD dust map 
and the Fermi 0.5 — 1 GeV IC template, the correlation 
coefficients are weighted by the mean of each template in 
the "bubble" region. The resulting model and the differ- 
ence maps with respect to the Fermi data, at different 
energy bands, are shown in Figure [131 The residuals are 
remarkably small. The spectrum is shown in Figure [Ml 
The Fermi 0.5 — 1 GeV IC template appears to contain 
a small fraction of tt^ gammas, but the spectral index 
is consistent with the predicted GALPROP IC component. 
The fitting coefficients and corresponding errors of each 
template are listed in Table [H 

By eye, the Fermi bubbles appear to possess north- 
south symmetry and are close to spatially uniform in 
intensity (with a hard cutoff at the bubble edges). To 
test these hypotheses more quantitatively, we split the 
Fermi bubble template into the inner bubble and outer 
shell templates {upper row of Figure [TS]) . or alternatively 
into the north and south bubble templates {lower row of 
Figure [15]). We then repeat the previous 5-template and 
4-template fitting procedure involving either the simple 
disk IC template or the Fermi 0.5 — 1 GeV IC template, 
but splitting the bubble template to either inner bub- 
ble and outer shell or north and south bubble templates. 
The goal is to identify variations in the intensity and 
spectral index between the bubble edge and interior, and 
the northern and southern bubbles. The resulting fitted 
spectra are shown in Figure [T6l And the corresponding 
fitting coefficients and errors of each template are listed 
in Table [5l Table [6l Table [71 and Table [HI respectively. 
In Figure [TTl we replace the simple disk model with the 
Haslam 408 MHz map as the IC template, and employ 
the SFD map, a uniform background, the Loop I tem- 
plate, and the double- lobed bubble {left panel) or the 
north and south bubble {right panel) templates in the fit- 
ting. Our conclusion is that the Fermi bubbles appear 
to be north-south symmetric and spatially and spectrally 
uniform, with a hard spectrum. This statement is largely 
independent of our choice of template for the disk IC 
emission. 



3.3. The ^Termi Bubbles'' vs the 'Termi Haze'' 
Dobler et al.l (|2010[ ) employed essentially the same 



template regression methods, claimed that the gamma- 



ray emission not accounted for by the known foregrounds 
could be well fitted by a bivariate Gaussian with (75 = 
25°, cFi = 15° . With the improved 1.6 yr data, the edges 
of the excess at high latitudes are seen to be quite sharp, 
and are not well-described by a Gaussian fall-off in in- 
tensity. However, the question of whether the excess is 
better modeled as an "egg" or a pair of "bubbles" is more 
subtle. 

The choice of t he bivariate Gaussian template by 
Dobler et al.l (|2010[ ) was intended to remove as much of 
the remaining gamma ray signal as possible, once the tt^ 
and soft IC emission had been regressed out, minimizing 
large-scale residuals. The fits were performed with only 
\h\ < 5 masked out. In this work, on the other hand, 
we have masked out all emission with \b\ < 30°; when 
attempting to subtract the disk-correlated emission, and 
delineating the "bubbles" template, our goal has been 
to isolate the sharp-edged features from the more slowly 
spatially varying emission, not to account for all the ob- 
served emission. 

This difference in approach can be seen in the non- 
negligible residuals around the inner Galaxy (|&| ^ 20°) 
in many of our maps (for example. Figure [13]). An at- 
tempt to fit all the residual emission simultaneously with 
a simple ten iplate may well requ ire a template closer to 
that used bv lDobler et al.l (|2010l ). rather than the bubble 
template. However, the sharp edges now visible in the 
data, and their alignment with the edges of the WMAP 
haze and ROSAT X-ray features (as we will discuss in 
Q, motivate us to consider the bubbles as originating 
from a distinct physical mechanism. While the bubbles 
probably do not consti t ute th e entire ''Fermi haze" dis- 
cussed bv lDobler et al.l (|2010[ ). they are certainly a major 
component of it, dominating the signal at high latitudes. 
The question of whether the remaining residual gamma 
rays represent a separate physical mechanism is an inter- 
esting one that we defer to future work. 

4. THE FERMI BUBBLES SEEN IN OTHER MAPS 

In this section, we compare the 1 — 5 GeV Fermi bubble 
with the ROSAT 1.5 keV soft X-ray map, the WMAP 
23 GHz microwave haze, and the Haslam 408 MHz map. 
The striking similarities of several morphological features 
in these maps strongly suggest a common physical origin 
for the Fermi bubbles, WMAP haze, and X-ray edges 
towards the GC (Figure [TEj). 

4.1. Comparison with ROSAT X-ray Features 

The ROSAT all-sky survey provides full-sky images 
with FWHM 12' at energies from 0.5-2 keV.^"^ We com- 
pare the morphology of the X-ray features in ROSAT l.b 
keV map with the edges of the Fermi bubbles in detail 
in Figure [191 The limb brightened X-ray features align 
with the edges of both the north and south Fermi bub- 
ble. Hints of the whole north bubble are also visible in 

^'^ http://hea-www.harvard.edu/rosat/rsdc.html 
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TABLE 1 

Coordinates defining the features shown in Figure [4] and Figure d 



Uniform 


SFD 


bubble 


north bubble 


south bubble 


inner bubble 


outer bubble 


disk 


loop I 


0.5- 


- 1.0 GeV - SFD 


1.0 


0.084 


1.0 


1.0 


1.0 


1.0 


1.0 


0.292 


1.0 




1.198 



TABLE 2 
Normalization factors for different templates. 



E range (GeV) 


Energy 


Uniform 


SFD dust 


simple disk 


whole bubble 


simple loop I 


0.3 


-0.5 


0.4 


1.681 ± 0.006 


1.201 ± 0.011 


0.689 ± 0.027 


0.035 ± 0.033 


0.487 ± 0.015 


0.5 


-0.9 


0.7 


1.365 ± 0.007 


1.279 ± 0.012 


0.608 ± 0.030 


0.211 ± 0.037 


0.475 ± 0.016 


0.9 


-1.7 


1.3 


1.141 ± 0.008 


1.179 ± 0.014 


0.503 ± 0.035 


0.321 ± 0.044 


0.405 ± 0.019 


1.7 


-3.0 


2.2 


1.034 ± 0.006 


0.876 ± 0.011 


0.393 ± 0.029 


0.436 ± 0.036 


0.376 ± 0.016 


3.0 


-5.3 


4.0 


0.881 ± 0.008 


0.554 ± 0.013 


0.420 ± 0.034 


0.353 ± 0.043 


0.249 ± 0.018 


5.3 


-9.5 


7.1 


0.731 ± 0.009 


0.322 ± 0.014 


0.282 ± 0.039 


0.343 ± 0.049 


0.208 ± 0.021 


9.5- 


-16.9 


12.7 


0.563 ± 0.010 


0.193 ± 0.015 


0.251 ± 0.044 


0.205 ± 0.055 


0.092 ± 0.023 


16.9- 


-30.0 


22.5 


0.507 ± 0.012 


0.128 ± 0.018 


0.191 ± 0.053 


0.263 ± 0.068 


0.125 ± 0.029 


30.0- 


-53.3 


40.0 


0.557 ± 0.015 


0.041 ± 0.021 


0.096 ± 0.064 


0.217 ± 0.083 


0.088 ± 0.036 


53.3- 


-94.9 


71.1 


0.628 ± 0.022 


0.020 ± 0.030 


0.183 ± 0.093 


0.251 ± 0.120 


-0.043 ± 0.048 


94.9- 


168.7 


126.5 


0.622 ± 0.030 


0.080 ± 0.043 


0.012 ± 0.125 


0.319 ± 0.162 


-0.091 ± 0.063 


168.7- 


300.0 


225.0 


0.436 ± 0.038 


0.174 ± 0.061 


-0.039 ± 0.158 


-0.015 ± 0.194 


0.083 ± 0.086 



TABLE 3 

Corresponding template fitting coefficients and errors in Figure [TJ] 



ROSAT, as well as two sharp edges in the south that 
trace the south Fermi bubble close to the disk. We show 
the ROSAT 1.5 keV map overplotted with the edges of 
the Fermi bubbles, the northern arc^ the ''donut" and 
the Loop I features in the right panels of Figure [191 The 
appearance of the X-ray edges in the ROSAT 1.5 keV 
map, coincident with the Fermi bubble edges, strongly 
supports the physical reality of these sharp edges. 
In Figure ^ we subtract the ROSAT 1.0 keV soft 



X-ray map from the 1.5 keV map to clean up the fore- 
ground. We find that the extended Loop I feature has a 
softer spectrum than the X-ray features associated with 
the bubble edges, and is largely removed in the differ- 
ence map {lower left panel of Figure [20]) . The residual 
features strikingly overlap with the edges of the Fermi 
bubbles {lower right panel). No other noticeable large 
scale features appear in the residual X-ray map which do 
not appear in the gamma-rays. 
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E range (GeV) 


Energy 


Uniform 


SFD dust 


whole bubble 0.5-1.0 GeV - SFD 


0.3 


-0.5 


0.4 


1.759 ± 0.006 


0.883 ± 0.012 


-0.026 ± 0.026 


1.181 ± 0.017 


0.5 


-0.9 


0.7 


1.446 ± 0.006 


0.905 ± 0.013 


0.008 ± 0.029 


1.275 ± 0.018 


0.9 


- 1.7 


1.3 


1.208 ± 0.008 


0.929 ± 0.016 


0.258 ± 0.034 


0.897 ± 0.020 


1.7 


-3.0 


2.2 


1.088 ± 0.006 


0.679 ± 0.012 


0.375 ± 0.029 


0.736 ± 0.017 


3.0 


-5.3 


4.0 


0.921 ± 0.007 


0.427 ± 0.014 


0.428 ± 0.034 


0.530 ± 0.019 


5.3 


-9.5 


7.1 


0.759 ± 0.008 


0.231 ± 0.015 


0.371 ± 0.039 


0.400 ± 0.021 


9.5- 


-16.9 


12.7 


0.580 ± 0.009 


0.131 ± 0.016 


0.269 ± 0.043 


0.270 ± 0.023 


16.9- 


-30.0 


22.5 


0.522 ± 0.011 


0.072 ± 0.016 


0.290 ± 0.053 


0.261 ± 0.020 


30.0- 


-53.3 


40.0 


0.565 ± 0.014 


0.015 ± 0.021 


0.225 ± 0.066 


0.141 ± 0.034 


53.3- 


-94.9 


71.1 


0.631 ± 0.021 


0.010 ± 0.033 


0.364 ± 0.096 


0.065 ± 0.053 


94.9- 


168.7 


126.5 


0.615 ± 0.029 


0.083 ± 0.048 


0.337 ± 0.127 


-0.032 ± 0.076 


168.7- 


300.0 


225.0 


0.441 ± 0.037 


0.149 ± 0.068 


-0.104 ± 0.137 


0.089 ± 0.097 



TABLE 4 
Corresponding template fitting coefficients and errors in Figure fT4] 



4.2. Comparison with WMAP Microwave Haze 

The WMAP haze is the residual remaining in WMAP 
microwave data after regressing out contributions from 
thermal dust, free- free, and "soft sy nchrotron" traced b y 
the Haslam 408 MHz radio survey (jHaslam et alJll982[ ). 
Therefore, it is by construction harder than the Haslam- 
correlated emission. We will show in this section that 
the WMAP synchrotron haze appears to be associated 
with the Fermi bubbles. 

Figure [2T] shows a detailed morphological comparison 
of the south Fermi bubble at 1 — 5 GeV with the southern 
part of the WMAP microwave haze at 23 GHz (K-band). 
The edge of the Fermi bubbles, marked in green dashed 
line in the top right and lower right panels, closely traces 
the edge of the WMAP haze. The smaller latitudinal 
extension of the WMAP haze may be due to the de- 
cay of the magnetic field strength with latitude. These 
striking morphological similarities between the WMAP 
microwave haze and Fermi gamma-ray bubble can be 
readily explained if the same electron CR population is 
responsible for both excesses, with the electron CRs in- 
teracting with the galactic magnetic field to produce syn- 
chrotron, and interacting with the ISRF to produce IC 
emission. 

In Figure [22l we show the difference maps between the 
23 GHz WMAP haze and the 33 GHz and 41 GHz haze 
maps. The difference maps contain no apparent features 
and indicate a common spectrum of different regions in- 
side the WMAP haze, suggesting a single physical origin 
for the bulk of the signal. We have reached the same 
conclusion for the Fermi bubbles in Figure [71 

Besides the similarity of the morphology, the relatively 
hard spectrum of the Fermi bubble also motivates a com- 
mon physical origin with the WMAP haze. We now 
provide a simple estimate of the microwave synchrotron 
and gamma-ray ICS signals from one single population of 
hard electrons distributed in the inner Galaxy, to demon- 
strate that the magnitudes and spectral indices of the two 
signals are consistent for reasonable parameter values. 

For a highly relativistic electron scattering on low 



energy phot ons, the spectrum of upsca tt ered photons 
is giv en by iBlumenthal fc Gouldl (|l97Qh (jCholis et al. 
200i, 



dN 



dEjdedt 4 



^^^^l!!^ (2g log g + (1 + 2q){l - q) 



eE^ 



+0.5(1 -g)(rg)V(l + rg))n(e), 
T = 4eEe/{meC^)\ q 



(2) 



^er(l-^^/^e)' 



e<E^ <£;er/(i + r). 

Here e is the initial photon energy, Ee is the electron en- 
ergy, E^ is the energy of the upscattered gamma ray, and 
n(e) describes the energy distribution of the soft photons 
per unit volume. Where F <C 1, in the Thomson limit, 
the average energy of the upscattered photons is given 

by, 

{E,} = (4/3)7^(e), (3) 

where 7 = Ee/nieC^ is the Lorentz boost of the electron. 
In the Klein-Nishina (KN) limit, F ^ 1, the spectrum 
instead peaks at the high energy end, and the upscattered 
photon carries away almost all the energy of the electron. 
Given a power law steady-state electron spectrum with 
spectral index 7, the spectral index of the IC scattered 
gamma rays is (7 + l)/2 in the Thomson limit, and 
7 + 1 in the extreme KN limit (e.g. IBlumenthal fc Gouldl 



I970I ). Photons in the Fermi energy range can be pro- 
duced by scattering of (9(10 — 100) GeV electrons on 
starlight, which is (marginally) in the KN regime, or by 
Thomson scattering of much higher-energy electrons on 
IR or CMB photons. Consequently, the spectral index 
of gamma rays might be expected to vary with latitude 
even if the electron spectral index is uniform, becoming 
harder at higher latitudes where scatterings in the Thom- 
son limit dominate. Closer to the disk, where much of 
the ISRF energy density is in starlight, IC scatterings are 
in neither the extreme KN nor Thomson limits, and the 
spectrum needs to be computed carefully. 

We consider a steady-state electron spectrum described 
by a power law, dN/dE oc E~^ , with energy cutoffs at 
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Template 



00.5-01.7GeV 



Residual map 




Fig. 11. — The models obtained from the multi-template fits, compared with the Fermi maps, in different energy bins. The left 
column shows the linear combination of the disk. Loop /, uniform, and bubble templates that provide the best fit to the Fermi maps after 
subtracting the best fit SFD dust template (shown in the middle column). The difference maps between the combined template and the 
data are shown in the right column. The template fitting is done for the region with |6| > 30° to avoid contaminations from the Galactic 
disk (shown with black dashed line in the right column residual maps). The subtraction of the model largely removes the features seen in 
the Ferm^i maps with |6| > 30°. We use the same gray scale for all the panels. We find that both the disk IC template and Loop I features 
fade off with increasing energy, but the bubble template is required for all the energy bands and does not fade off with increasing energy. 
The oversubtraction in the residual maps, especially at lower energy bins, is due to the simple disk IC model, which is not a good template 
across the entire disk. However, outside the masked region, the residual maps are consistent with Poisson noise without obvious large scale 
features. 
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Fig. 12. — Correlation spectra for the 5-template fit employing a simple disk model for the IC (and to a lesser degree bremsstrahlung) 
emission from supernova-shock-accelerated electrons (see ^3.2.2|) . The SFD-correlated spectrum is shown by the red short-dashed line which 
roughly traces tt*-* emission (the gray dashed line indicates a GALPROP prediction for tt*-* emission). The disk-correlated emission is shown 
by the green dashed line, which traces the soft IC (gray triple-dot-dashed line) and bremsstrahlung (gray dot-dashed line) component. 
The spectrum of the uniform emission, which traces the isotropic background (including possible cosmic-ray contamination), is shown as 
a dotted brown line. The solid orange line indicates the spectrum of emission correlated with Loo-p /, which has a similar spectrum to the 
disk-correlated emission. Finally, the blue dot-dashed line shows the spectrum correlated with the Fermi bubble template. The bubble 
component has a notably harder (consistent with fiat) spectrum than the other template-correlated spectra, and the models for the various 
emission mechanism generated from GALPROP, indicating that the Fermi bubbles constitute a distinct component with a hard spectrum. 
The fitting is done over the |6| > 30° region. Note that these GALPROP "predictions" are intended only to indicate the expected spectral 
shape for these emission components, for reference. The correlation coefiicients for the SFD map and simple disk model are multiplied by 
the average value of these maps in the bubble region (defined by the bottom right panel of Figure [S] with a |6| > 30° cut) to obtain the 
associated gamma-ray emission; see ^3.2.2l for details, and Table [2] for a summary of the normalization factors. 

the corresponding synchrotron and IC spectra. The 
WMAP haze was esti mated to have a spectrum ly oc 



0.1 GeV and 1000 GeV. The choice of high-energy cut- 
off is motivated by the local mea surement of the c osmic 
ray electron spectrum by Fermi (JAbdo et al.luOOQf ). We 
consider a region ^ 4 kpc above the Galactic center, as 
an example (and since both the WMAP haze and Fermi 
bubbles are reasonably well measured there), and em- 
ploy the model for the ISR F used in GALPROP version 
50p (jPorter fc StrondboOSh at 4 kpc above the GC. We 
normalize the synchrotron to the approxim ate value mea- 



-/5 



, P = 0.39 - 0.67 (JDobler fc Finkbeinerll2008[ ). cor 



sured by WMAP in the 23 GHz K-band (jHooper et al 
20071 ). ~ 25° below the Galactic plane, and compute 



responding approximately to an electron spectral index 
of 7 ^ 1.8 — 2.4; Figure [23] shows our results for a mag- 
netic field of 10 /iG and 5 /iG at 4 kpc above the GC, 
and electron spectral indices 7 = 1.8 — 3. We find good 
agreement in the case of a ~ 2 — 2.5, consistent with the 
spectrum of the WMAP haze. 
In the default GALPROP exponential model for the 



Galactic magnetic field, \B\ 



\Bo\e- 



-Z I Zs 



with scale 
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Fig. 13. — The same as Figure [TT] but using the Fermi 0.5 — 1 GeV residual map after subtracting emission correlated with the SFD 
dust map (to remove the tt*-" gammas) as a template for low-latitude IC emission (originating primarily from scatterings on starlight, and 
likely involving softer supernova-shock- accelerated electrons). The left column shows the best fit combined template including the residual 
0.5 — 1 GeV map, uniform, SFD dust, and bubble templates. The middle column shows the Fermi-LAT data in different energy bins. The 
difference maps between the combined template and the data are shown in the right column. The template fitting is done over the region 
with |6| > 30° (shown with black dashed line in the right column residual maps). We use the same grey scale for all the panels. We find 
that the bubble template does not fade away with increasing energy. The 4-template fit works extremely well, at |6| > 30° the residual 
maps are consistent with Poisson noise without obvious large scale features, and there are no obvious sharp features in the data closer to 
the disk. 
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Fig. 14. — Same as Figure fT2l but correlation spectra for the 4-template fit employing the Fermi 0.5 — 1 GeV residual map (after 
subtracting the SFD dust) as a template for the starlight IC. The line style is the same as Figure fT2] Again, we find that the spectrum 
correlated with the Fermi bubble template (blue dot-dashed line) is harder (consistent with fiat in E'^dN/dE) than the spectra correlated 
with the other templates, and the models for the various emission mechanism generated from GALPROP, indicating that the Fermi bubbles 
constitute a distinct gamma-ray component with a hard spectrum. The fitting is done for |6| > 30°. As in Figure fT2l the correlation spectra 
have been normalized to a reference region; see ^3.2.2l for details, and Table [2] for a summary of the normalization factors. 



height Zs ^ 2 kpc, this field strength would correspond to 
5o ^ 30 — 40 /iG or even higher. This value is consider - 
ably larger than commonly used (e.g. iPage et al J 120071 ). 
However, models with a n on-exponential halo n i agneti c 
field, as d i scusse d by e.g. lAlvarez-Muhiz et al.l (J2002|); 
Sun et all (|2008h . can have - 10 /iG fields weh off the 
plane. 

We note also that the extrapolated value of Bq re- 
quired to obtain good agreement between the IC and syn- 
chrotron amplitudes, in t he exponential mode l, is some- 
what higher than found bv lPobler et al.l (|2010[ ). who per- 
formed the comparison at 2 kpc. This apparent discrep- 
ancy originates from the fact that in the haze latitudi- 
nal profile given by iHooper et al.l (|2007[ ). the emission 
falls off rapidly with latitude for > h > —15°, but 
then plateaus at 6 ~ —15 — 35°, contrary to expectations 



based on a B-field profile exponentially falling away from 
z = 0. This suggests either that the magnetic field inside 
the bubble does not fall exponentially with \z\ inside the 
bubbles, or that the WMAP haze contains a significant 
free- free component at high latitude. 

4.3. Evidence of a ^ 700 GeV Electron Excess? 

In Figure [24l we calculate the gamma-ray spectrum 
from IC scattering, using the standard ISRF model taken 
from GALPROP - as in Figure [23l but with a different en- 
ergy range for the electron CRs. An electron CR popula- 
tion with a hard low-energy cutoff at about 500 GeV can 
fit the Fermi bubble spectrum better than a single power 
law extending from 0.1 — 1000 GeV, due to the downturn 
in the spectrum in the lowest energy bin. Even a rather 
hard {dN/dE ^ E~'^) power law component at 300 — 500 
GeV produces a long tail at low energies. Interestingly, 
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Fig. 15. — Top row: We split our Fermi bubble template (shown in the bottom right panel of Figure [3]) into two components for template 
fitting: an interior template (top left) and a shell template (top right) with uniform intensity, in order to reveal any potential spectrum 
difference with the template fitting technique. Bottom, row: We split the Fermi bubble template into north bubble template (bottom left) 
and south bubble template (bottom right). If the two bubbles have the same origin, they should not only have similar morphologies but 
also consistent spectra. 



this preferred 500 — 700 GeV energy range is rather close 
to the peak in local e^ - f e~ cosmic rays observed by 
ATIC (jChang et alJboOSl )^^. We note that although the 
estimated error bars of energy lower than 1 GeV mildly 
depend on the templates we use in the fitting procedure, 
the fall-off of the bubble intensity in the lowest energy 
bins is robust. Figure [23] shows the same analysis for 
500-900 GeV electrons with two different templates for 
the disk IC emission, for the bubble interior and shell 
separately, and for the north and south bubbles. In all 
cases the same cut-off in the spectrum below 1 GeV is 
observed. 

In this case, while the lowest energy bin in the gamma- 
rays is better fitted, the lack of low-energy CRs means 

^^ However, note that the large peak observed by ATIC -2 and 
ATIC -4, appears to be in confiict with the Ferm^i measurement of 
the spectrum, which contains no such feature. 



that synchrotron can contribute to the WMAP haze only 
at the sub-percent level, unless the magnetic field in the 
inner galaxy is extreme (and even then the spectrum 
does not reproduce the observations). In this case an 
alternate explanation for the WMAP haze would need 
to be c onsidered. As suggested initially in iFinkbeiner 
( 2004al ). the WMAP haze could originate from free- free 
emission (thermal bremsstrahlung), and this would ex- 
plain the lack of a clear haze signal in WMAP polariza- 
tion maps. However, the spectrum of the haze is some- 
what softer than generally expected from free-free emis- 
sion; the gas temperature required is also thermally un- 
stable, requiring a significant energy injection (^ 2q54-55 



ergs) to maintain its tempera ture (jHooper et al.l 12007 : 
McQuinn fc Zaldarriagall2010h . 



4.4. Gamma-ray power and e cosmic ray density 
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E range (GeV) 


Energy 


Uniform 


SFD dust 


simple disk 


inner bubble 


outer bubble 


simple loop I 


0.3 


-0.5 


0.4 


1.681 ± 0.006 


1.201 ± 0.011 


0.683 ± 0.027 


0.071 ± 0.042 


0.004 ± 0.040 


0.490 ± 0.015 


0.5 


-0.9 


0.7 


1.365 ± 0.007 


1.279 ± 0.012 


0.607 ± 0.030 


0.215 ± 0.048 


0.207 ± 0.045 


0.475 ± 0.016 


0.9 


- 1.7 


1.3 


1.142 ± 0.008 


1.179 ± 0.014 


0.498 ± 0.036 


0.354 ± 0.057 


0.293 ± 0.053 


0.407 ± 0.019 


1.7 


-3.0 


2.2 


1.034 ± 0.006 


0.876 ± 0.011 


0.403 ± 0.029 


0.370 ± 0.047 


0.491 ± 0.045 


0.373 ± 0.016 


3.0 


-5.3 


4.0 


0.880 ± 0.008 


0.554 ± 0.013 


0.426 ± 0.035 


0.307 ± 0.056 


0.393 ± 0.054 


0.247 ± 0.018 


5.3 


-9.5 


7.1 


0.731 ± 0.009 


0.322 ± 0.014 


0.284 ± 0.039 


0.330 ± 0.064 


0.354 ± 0.061 


0.207 ± 0.021 


9.5- 


-16.9 


12.7 


0.562 ± 0.010 


0.193 ± 0.015 


0.265 ± 0.044 


0.100 ± 0.070 


0.293 ± 0.069 


0.087 ± 0.023 


16.9- 


-30.0 


22.5 


0.506 ± 0.012 


0.128 ± 0.018 


0.201 ± 0.054 


0.182 ± 0.088 


0.335 ± 0.087 


0.121 ± 0.029 


30.0- 


-53.3 


40.0 


0.556 ± 0.015 


0.041 ± 0.021 


0.098 ± 0.065 


0.190 ± 0.106 


0.243 ± 0.107 


0.087 ± 0.036 


53.3- 


-94.9 


71.1 


0.627 ± 0.022 


0.020 ± 0.030 


0.187 ± 0.093 


0.206 ± 0.151 


0.294 ± 0.153 


-0.045 ± 0.048 


94.9- 


168.7 


126.5 


0.620 ± 0.030 


0.081 ± 0.043 


0.037 ± 0.129 


0.180 ± 0.209 


0.431 ± 0.207 


-0.098 ± 0.064 


168.7- 


300.0 


225.0 


0.435 ± 0.038 


0.179 ± 0.062 


-0.065 ± 0.155 


0.145 ± 0.242 


-0.178 ± 0.206 


0.097 ± 0.086 



TABLE 5 
Corresponding template fitting coefficients and errors in upper left panel of Figure [16] 



E range (GeV) 


Energy 


Uniform 


SFD dust 


simple disk 


north bubble 


south bubble 


simple loop I 


0.3 


-0.5 


0.4 


1.679 ± 0.006 


1.205 ± 0.011 


0.683 ± 0.027 


-0.042 ± 0.050 


0.069 ± 0.037 


0.496 ± 0.015 


0.5 


-0.9 


0.7 


1.363 ± 0.007 


1.282 ± 0.012 


0.603 ± 0.030 


0.140 ± 0.056 


0.241 ± 0.041 


0.483 ± 0.017 


0.9 


-1.7 


1.3 


1.138 ± 0.008 


1.187 ± 0.015 


0.493 ± 0.035 


0.171 ± 0.065 


0.385 ± 0.049 


0.421 ± 0.020 


1.7 


-3.0 


2.2 


1.030 ± 0.007 


0.883 ± 0.012 


0.384 ± 0.029 


0.288 ± 0.054 


0.503 ± 0.041 


0.391 ± 0.016 


3.0 


-5.3 


4.0 


0.881 ± 0.008 


0.553 ± 0.013 


0.421 ± 0.035 


0.368 ± 0.064 


0.346 ± 0.049 


0.247 ± 0.019 


5.3 


-9.5 


7.1 


0.731 ± 0.009 


0.321 ± 0.014 


0.283 ± 0.039 


0.354 ± 0.072 


0.337 ± 0.055 


0.207 ± 0.022 


9.5- 


-16.9 


12.7 


0.562 ± 0.010 


0.195 ± 0.016 


0.248 ± 0.044 


0.157 ± 0.078 


0.230 ± 0.063 


0.097 ± 0.024 


16.9- 


-30.0 


22.5 


0.509 ± 0.012 


0.124 ± 0.018 


0.196 ± 0.053 


0.353 ± 0.100 


0.216 ± 0.076 


0.116 ± 0.030 


30.0- 


-53.3 


40.0 


0.556 ± 0.015 


0.042 ± 0.021 


0.094 ± 0.064 


0.190 ± 0.117 


0.232 ± 0.097 


0.092 ± 0.038 


53.3- 


-94.9 


71.1 


0.628 ± 0.022 


0.020 ± 0.030 


0.184 ± 0.093 


0.264 ± 0.164 


0.243 ± 0.140 


-0.045 ± 0.049 


94.9- 


168.7 


126.5 


0.614 ± 0.030 


0.091 ± 0.044 


0.004 ± 0.126 


0.110 ± 0.208 


0.450 ± 0.195 


-0.068 ± 0.067 


168.7- 


300.0 


225.0 


0.429 ± 0.039 


0.186 ± 0.064 


-0.058 ± 0.160 


-0.169 ± 0.254 


0.061 ± 0.222 


0.107 ± 0.092 



TABLE 6 
Corresponding template fitting coefficients and errors in upper right panel of Figure [161 



E range (GeV) 


Energy 


Uniform 


SFD dust 


inner bubble 


outer bubble 0.5 - 1.0 GeV - SFD 


0.3 


-0.5 


0.4 


1.759 ± 0.006 


0.883 ± 0.012 


0.011 ± 0.035 


-0.069 ± 0.037 


1.180 ± 0.017 


0.5 


-0.9 


0.7 


1.446 ± 0.006 


0.905 ± 0.013 


0.005 ± 0.039 


0.012 ± 0.041 


1.275 ± 0.018 


0.9 


-1.7 


1.3 


1.208 ± 0.008 


0.929 ± 0.016 


0.284 ± 0.047 


0.230 ± 0.048 


0.896 ± 0.020 


1.7 


-3.0 


2.2 


1.088 ± 0.006 


0.680 ± 0.012 


0.309 ± 0.039 


0.444 ± 0.041 


0.737 ± 0.017 


3.0 


-5.3 


4.0 


0.921 ± 0.007 


0.427 ± 0.014 


0.399 ± 0.047 


0.459 ± 0.049 


0.530 ± 0.019 


5.3 


-9.5 


7.1 


0.759 ± 0.008 


0.231 ± 0.015 


0.362 ± 0.053 


0.382 ± 0.056 


0.400 ± 0.021 


9.5- 


-16.9 


12.7 


0.580 ± 0.009 


0.131 ± 0.016 


0.195 ± 0.057 


0.351 ± 0.063 


0.271 ± 0.023 


16.9- 


-30.0 


22.5 


0.523 ± 0.011 


0.069 ± 0.015 


0.218 ± 0.071 


0.368 ± 0.079 


0.261 ± 0.018 


30.0- 


-53.3 


40.0 


0.565 ± 0.012 


0.014 ± 0.015 


0.202 ± 0.088 


0.249 ± 0.096 


0.144 ± 0.025 


53.3- 


-94.9 


71.1 


0.631 ± 0.021 


0.010 ± 0.033 


0.337 ± 0.128 


0.394 ± 0.140 


0.065 ± 0.053 


94.9- 


168.7 


126.5 


0.614 ± 0.029 


0.084 ± 0.048 


0.232 ± 0.164 


0.455 ± 0.190 


-0.029 ± 0.076 


168.7- 


300.0 


225.0 


0.440 ± 0.037 


0.152 ± 0.069 


0.021 ± 0.195 


-0.262 ± 0.175 


0.088 ± 0.096 



TABLE 7 
Corresponding template fitting coefficients and errors in lower left panel of Figure [16] 



In order to estimate the total gamma-ray power emit- 
ted by the bubbles, we must estimate the surface 
brightness integrated over energy, the solid angle sub- 
tended, and the distance (suitably averaged). From 
Figure [12] we take the intensity to be E'^dN/dE = 
3 X IQ-'^GeVcm-^s-^sr-^ from 1 - 100 GeV, integrat- 
ing to 1.4x10"^ GeV cm~^s~^sr~^. The bubble template 
used in our analysis (Figure [3|) subtends 0.808 sr, yield- 
ing a total bubble flux of 1.13x10"^ GeVcm-^s-^ To 
obtain an average distance for the emission, we approx- 
imate the bubbles as 2 spheres centered at 6 = ±28°, 



and directly above and below the Galactic center. For 
a Sun-GC distance of 8.5 kpc, this implies a distance of 
9.6 kpc, and a total power (both bubbles) in the 1 — 100 
GeV band of 2.5 xlO'^^ GeV/s or 4.0 x 10^^ erg/s, which 
is ~ 5% of the total Galactic gamma-ray luminosity be- 






tween 0.1 - 100 GeV (jStrong et al.ll201C 

The electron cosmic-ray density in the bubbles required 
to generate the observed gamma rays, at any given en- 
ergy, depends strongly on the assumed electron spec- 
trum. However, typically the required values are com- 
parable to the locally measured electron CR density. For 
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E range (GeV) 


Energy 


Uniform 


SFD dust 


north bubble 


south bubble 0.5 - 1.0 GeV - SFD 


0.3 


-0.5 


0.4 


1.759 ± 0.006 


0.882 ± 0.012 


-0.011 ± 0.047 


-0.031 ± 0.029 


1.180 ± 0.017 


0.5 


-0.9 


0.7 


1.446 ± 0.006 


0.905 ± 0.013 


-0.009 ± 0.052 


0.014 ± 0.032 


1.276 ± 0.018 


0.9 


- 1.7 


1.3 


1.207 ± 0.008 


0.931 ± 0.016 


0.201 ± 0.061 


0.277 ± 0.038 


0.899 ± 0.020 


1.7 


-3.0 


2.2 


1.087 ± 0.006 


0.681 ± 0.013 


0.331 ± 0.050 


0.391 ± 0.033 


0.738 ± 0.017 


3.0 


-5.3 


4.0 


0.922 ± 0.007 


0.424 ± 0.014 


0.500 ± 0.060 


0.399 ± 0.039 


0.527 ± 0.019 


5.3 


-9.5 


7.1 


0.760 ± 0.008 


0.228 ± 0.015 


0.440 ± 0.067 


0.342 ± 0.045 


0.397 ± 0.021 


9.5- 


-16.9 


12.7 


0.580 ± 0.009 


0.132 ± 0.016 


0.231 ± 0.072 


0.286 ± 0.051 


0.272 ± 0.023 


16.9- 


-30.0 


22.5 


0.524 ± 0.011 


0.068 ± 0.015 


0.406 ± 0.091 


0.236 ± 0.060 


0.257 ± 0.008 


30.0- 


-53.3 


40.0 


0.565 ± 0.014 


0.014 ± 0.015 


0.227 ± 0.105 


0.222 ± 0.077 


0.143 ± 0.015 


53.3- 


-94.9 


71.1 


0.631 ± 0.021 


0.011 ± 0.033 


0.339 ± 0.151 


0.376 ± 0.115 


0.067 ± 0.054 


94.9- 


168.7 


126.5 


0.608 ± 0.029 


0.091 ± 0.049 


0.067 ± 0.189 


0.470 ± 0.162 


-0.013 ± 0.077 


168.7- 


300.0 


225.0 


0.437 ± 0.038 


0.156 ± 0.070 


-0.216 ± 0.229 


-0.066 ± 0.160 


0.095 ± 0.098 



TABLE 8 
Corresponding template fitting coefficients and errors in lower right panel of Figure fT6] 



example, for the model in the first panel of Figure [24] 
{dN/dE oc E-^ for 500 GeV < £; < 700 GeV), the in- 
ferred bubble electron density is ~ 10 x greater than the 
local electron density (as measured by Fermi) at an en- 
ergy of 500 GeV. For a representative model from the 
first panel of Figure [23] {dN/dE ex E'^-^ for 0.1 GeV 
< E < 1000 GeV, with a 10 fiG magnetic field generat- 
ing the WMAP Haze via synchrotron), at 500 GeV the 
bubble electron density is a factor ~ 2 x greater than the 
local density. 

5. INTERPRETATION 

As discussed in JDobler et al.l (|2010[ ). the Fermi bubbles 
seem most likely to originate from IC scattering, since 
the required electron CR population can also naturally 
generate the WMAP haze as a synchrotron signal. The 
ROSAT X-ray measurements suggest that the bubbles 
are hot and hence underdense regions, and thus argue 
against the gamma rays originating from bremsstrahlung 
or TT^ decay. 

Even though the material in the bubbles is likely 
high pressure, it is also probably very hot (~ 10^ K) 
and has lower gas density than the ambient ISM. This 
would explain why the ROSAT 1.5 keV map shows a 
"cavity" of soft X-rays toward the center of the Fermi 
bub ble structure, like the X-ray cavity in galaxy clus- 
ters (jMcNamara fc Nulsenll2007[ ). especially in the north 
Fermi bubble. Furthermore, allowing the Fermi bubbles 
to have lower density than the ambient medium means 
they would experience a buoyant force moving the bub- 
ble material away from the GC (see §6.2.3l for further dis- 
cussion) , which may help generate the observed morphol- 
ogy. Because tt^ and bremsstrahlung gamma-ray emis- 
sion both scale as the CR density x the gas density, an 
underdense region cannot be brighter unless the cosmic 
ray densities are greatly increased to compensate; for 
protons, in particular, the propagation lengths are great 
enough that a proton overdensity cannot reasonably ex- 
plain the sharp bubble edges observed in the data, if the 
bubbles are in a steady state. 

If the bubbles are expanding rapidly and highly accel- 
erated protons responsible for the gamma-ray emission 



are trapped behind shock fronts, then sharp edges for 
the Fermi bubbles could occur naturally. However, in 
the presence of such a shock, electrons would also be 
accelerated, and would generally produce more gamma- 
rays than the protons via ICS (since the cooling time for 
electron CRs is much shorter than the cooling time for 
proton CRs of comparable energy). 

It might be thought that the presence of a bright X- 
ray edge could lead to a sharp edge in the gamma-ray 
signal, via IC scattering of electron CRs on the X-ray 
photons. In the Thomson limit, the energy of IC scat- 
tered photons is of order {Te/'2)Ee^ with Fg = AEeE^/ml 
(where E^ and E^ are the initial photon and electron en- 
ergies, respectively, and rUe is the electron mass), and 
the scattering cross section is independent of the ini- 
tial electron and photon energies. Thus a higher-energy 
photon population, leading to a larger value of Fg, al- 
lows IC gamma-rays at a given energy to originate from 
lower-energy electrons, which are much more abundant 
for typical electron spectra with dN/dE ~ E~^ ^ 7^2. 
However, in the Klein-Nishina regime where Fg ^ 1 this 
picture changes: the energy of scattered photons is de- 
termined mostly by the energy of the initial electron, 
and the cross section scales as l/Fg. Scatterings of ~ 50 
GeV electrons already produce ~ 10 GeV gamma rays; 
when compared to the scattering of 10 GeV electrons on 
X-rays in the extreme KN limit, the KN suppression of 
the cross section in the latter case more than counteracts 
the greater abundance of ~ 10 GeV electrons, unless the 
electron spectrum is very soft (which is inconsistent with 
the observed signal). 

Thus hard UV or X-rays in the bubbles, which might 
be naturally expected in a high-temperature region, 
would not make the IC spectrum harder, and IC scat- 
tering on these photons is subdominant to IC scattering 
on the usual ISRF for the electron energies in question, 
unless the X-ray photon number density is much greater 
than the starlight photon number density. Furthermore, 
there is no reason to think the bubbles contains more 
< 1 eV photons, at least not in a region with a well de- 
fined spatial edge. Thus the sharp edges of the Fermi 
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Fig. 16. — Same as Figure [12] and Figure [M] but splitting the Fermi bubble template into two components for template fitting. The line 
styles are the same as Figure fT2] Top row: Using the simple disk model as the IC template. In the left panel, we split the previous bubble 
template into bubble interior and bubble shell templates (see Figure fTsl for the templates). The correlation coefficients of the 6-template 
fit involving the two bubble templates are shown. The purple dash-dotted line and blue triple-dot-dashed line are for the inner bubble and 
the outer shell template respectively. The two templates have a consistent spectrum which is significantly harder than the other templates, 
indicating the bubble interior and the bubble shell have the same distinct physical origin. In the right panel, we split the bubble template 
into north and south bubbles. As we include the Loop I template (which has a softer spectrum) in the north sky for regression fitting, 
the north bubble has a slightly harder spectrum than the south bubble. Again, both of the templates have harder spectra than any other 
components in the fit. Bottom row: Employing the Fermi 0.5 — 1 GeV residual map (after subtracting the SFD dust) as a template for 
the starlight IC In the left panel, we split the bubble template into bubble interior and bubble shell templates. In the right panel, we 
split the bubble template into north bubble and south bubble templates. As in Figure [T2] the correlation spectra have been normalized to 
a reference region; see ^3.2.2l for details, and Table [2] for a summary of the normalization factors. 
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Fig. 17. — Same as Figure fT2] but using the Haslam 408 MHz map instead of the simple disk model as the IC template. The Haslam 
map contains a bright feature associated with Loop I (see Figure [TS]) and is dominated by synchrotron emission from softer electron CRs, 
of energies around 1 GeV; it is not an ideal tracer of IC emission which depends on both electron and ISRF distribution. The resulting 
best-fit spectrum for the bubble template remains harder than the other components, but with enhanced lower energy (< 2 GeV) correlation 
coefficients compared to Figure [12] and Figure [M] The right panel is the same as the left panel, but with the bubble template divided 
into north and south bubbles (see ^3.2.2l for more discussion). As in Figure [T2] the correlation spectra have been normalized to a reference 
region; see ^3.2.2l for details, and Table [2] for a summary of the normalization factors. 



bubbles, and the non- uniformity in the emissivity, most 
Hkely arise from the electron CR density rather than the 
photon density. The presence of similar sharp edges in 
the WMAP haze (at \h\ < 30°) supports this hypothesis, 
if the WMAP haze is attributed to synchrotron radiation 
from the electron CRs. 

Similarly, the elongated shape of the Fermi bubble 
structures perpendicular to the Galactic plane suggests 
that the electron CR distribution itself is extended per- 
pendicular to the plane. The Fermi bubble morphol- 
ogy is a strong argument against the possibility that the 
WMAP haze originates from a disk-like electron distri- 
bution with significant lon ^tudinal variation of the mag- 
netic field, as suggested bv iKaplinghat et al.l (J2QQ9i. 

The limb brightening of the X-rays in the ROSAT dsidi 
(as shown in Figure [19]), and the flat intensity profile 
of the Fermi bubbles, suggest the presence of a shell 
or shock, with increased electron CR density, coincid- 
ing with a hot thermal plasma. If the ambient medium 
several kpc above and below the GC were neutral, then 
bubbles of ionized gas could produce a void in the Hi 
map (Figure [26]). We see no evidence for features aligned 
with the bubbles in these maps, suggesting that the H i 
map in this part of the sky is dominated by disk emis- 
sion, and has nothing to do with the bubbles. If the 
bubbles are in a static state, the bubble edges should 
have lower temperature than the bubble interior and thus 
higher gas density, although shocks or MHD turbulence 
might lead to higher temperatures at the bubble wall. 



The X-rays in ROSAT vad^y be thermal bremsstrahlung 
emission, so the emissivity is proportional to the ther- 
mal electron density x ion density. They could also 
arise from charge exchange reactions occurring when the 
high-speed gas in the bub bles collides w ith the denser 
gas at the bubble edge (see ISnowdenll2QQ9l . and reference 
therein); this mechanism could explain the pronounced 
limb brightening of the X-rays. As an alternative expla- 
nation, the ROSAT^i-idiy feature might be synchrotron 
emission from very high energy electron CRs. Typically, 
though, one needs ^ 50 TeV (~ 5 TeV) electrons with 
~ 10 fiG (1 mG) magnetic field to produce ~ 1 keV 
synchrotron photons. 

The Fermi bubble features do not appear to be as- 
sociated wit h Loop /, a giant radio loop spanning over 
100 degrees (JLarge et al.lll962[ ). which is thought to be 
generated from the local Sco-Cen OB association. De- 
tections of Loop I in h i gh-en ergy gamma-rays have been 



claimed bvlBhat et al.l (Il985l) and also recently by Fermi 
( Casandjian et al.ll2009[ ): we have also discussed its pres- 



ence in this work (see §3.2.2p . 

The Loop I gamma-rays may be the IC counterpart 
of the synchrotron emission seen in the Haslam 408 MHz 
map, although some of the emission might be tt^ gammas 
associated with Hi (Figure [26]). We compare structures 
identified from the Fermi 1 — 5 GeV maps with Loop I 
features in the Haslam 408 MHz map in the top row of 
Figure [iHl and see that the Fermi bubbles are spatially 
distinct from the arcs associated with Loop /; as we have 
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Fig. 18. — Comparison of the Fermi bubbles with features in other maps. Top left: Point-source subtracted 1 — 5 GeV Fermi-LAT 1.6 
yr map, same as the lower left panel of Figure [3] with north and south bubble edges marked with green dashed line, and north arc in blue 
dashed line. The approximate edge of the Loop I feature is plotted in red dotted line, and the "donut" in purple dot-dashed line. Top 
right: The Haslam 408 MHz map overplotted with the same red dotted line as the top left panel. The red dotted line remarkably traces 
the edge of the bright Loop I feature in the Haslam soft synchrotron map. Bottom left: the ROSAT 1.5 keV X-ray map is shown together 
with the same color lines marking the prominent Fermi bubble features. Bottom right: WMAP haze at K-band 23 GHz overplotted with 
Fermi bubble edges. The i^OiS'AT X-ray features and the WMAP haze trace the Fermi bubbles well, suggesting a common origin for these 
features. 

shown in Figure [161 the Loop I correlated emission also 
has a softer spectrum than the Fermi bubble emission. 
The Loop I feature in the ROSAT vnd^^p similarly has a 
softer spectrum than the limb-brightened X-ray bubble 
edges: as shown in Figure [20l when a low-energy map 
is subtracted from a higher-energy map in such a way 
that Loop I vanishes, the bubble edges remain bright. 
We also see additional shell structures which follow the 
Fermi bubble edges and the northern arc in the Haslam 
408 MHz map {top row of Figure [26]) . 

The Fermi bubbles are morphologically and spectrally 
distinct from both the tt^ emission and the IC and 
bremsstrahlung emission from the disk electrons. As we 
have shown in Figure [12] to Figure [TT] the Fermi bub- 
bles have a distinctly hard spectrum, dN^/dE ^ £^~^, 



with no evidence of spatial variation across the bub- 
bles. As shown in Figure [23] an electron population 
with dNe/dE ~ ^-2-2.5 ^g req^jire^j to produce these 
gamma rays by IC scattering: this is comparable to the 
spectrum of electrons a ccelerated by supernov a shocks or 
polar cap acceleration (JBiermann et al.ll2QlQ[ ). However, 
diffusive propagation and cooling would be expected to 
soften the spectrum, making it difficult to explain the 
Fermi bubbles by IC scattering from a steady-state pop- 
ulation of these electrons (a single brief injection of elec- 
trons with dN/dE ^ E~'^ could generate a sufficiently 
hard spectrum for the bubbles if there was a mechanism 
to transport them throughout the bubble without sig- 
nificant cooling). The facts strongly suggest that a dis- 
tinct electron component with a harder spectrum than 
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Fig. 19. — Top left: The i^O.S'AT X-ray haze features compared with the Fermi bubbles' morphology. Top row: The X-ray features of 
ROSAT hand 6 and 7 in the north sky towards to the GC (left panel) compared with Fermi north bubble overplotted with green dashed 
line, northern arc feature in blue dashed line, and Loop I feature in red dotted line (right panel). Bottom row: Same as top row but for 
the south bubble features. 

steady-state SNR generated electrons is responsible for 
the Fermi bubbles and associated signals in the WMAP 
and ROSAT dsitsi. 

It has been suggested that a large population of faint 
millisecond pulsars (MSPs) in the Milky Way halo could 
contribute to the WMAP and Fermi haze signals, via 
both pulsed ^amma-r ay emission and e+e~ production 
( Malvshevet alJboiol ). With a halo population of 3 x 10"^ 
MSPs, roughly half of the spin-down power going into 



disk population of pulsars, either peaking in the GC or 
peaking at small Galactocentric radius but goin^ to zero 
in the GC (jKaplinghat et alfcOOQclHarding fc Abazajian 
2Q1Q| ): such models have difficulty explaining the spheri- 



pairs and ~ 10% going to pulsed gamma-rays, con- 
sistency with both the WMAP data and the ffist-year 
Fermi photon data is possible; however, this model does 
not immediately explain either the rather sharp edge in 
the distribution of gamma-ray emission, or the features 
in the ROSAT X-ray data (the same can of course be 
said for models which generate the hazes via dark mat- 
ter annihilation or decay). Other attempts to explain 
the WMAP haze with pulsars have generally employed a 



cal morphology of the haze. Furth ermore, as pointed out 
by iMcQuinn fc Zald arriagal (120101 ) , the regression of the 
408 MHz Haslam map should remove much of the con- 
tribution from young pulsars in the WMAP data, since 
young pulsars have a similar spatial distribution to su- 
per novae. 

How could the electron CRs possess the same hard 
spectrum everywhere within the bubble and extend up 
to 10 kpc, while experiencing a steep fall-off at the bub- 
ble edge? A large population of CRs might be entrained 
in large scale Galactic outflows from the GC and enrich 
the bubbles (see more discussion in §6.3p . CRs could be 
produced along with jets, or shock accelerated CRs from 
magnetic reconnection inside the bubble or near its sur- 
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Fig. 20. — Top left: The i?0*S'AT X-ray features in Band 6 and 7. Top right: The same region as the left panel, but for the ROSAT X-ray 
map in band 5. Bottom left: The residual X-ray features after subtracting the top right softer Band 5 map from the top left harder Band 
6 and 7 map. Bottom right: The same as the bottom, left panel, overplotted with the Fermi bubble features, the northern arc, and Loop I 
features. The residual X-ray features with harder spectrum than the diffuse Loop I feature align well with the Fermi bubble structures. 



face (see more discussion in ^. However, it is challeng- 
ing to produce a fiat intensity profile for the bubble inte- 
rior with a sharp edge. The ambient gas should be com- 
pressed to a higher density on the shell by shocks (proba- 
bly also enhancing the magnetic field), and brighter syn- 
chrotron emission on the shell would then be expected, 
but the haze emission observed in WMAP is not limb- 
brightened and shows no evidence for a shell of finite 
thickness (although we do see shell structure in the X- 
rays). A cartoon picture summarizing the morphology of 
the Fermi bubbles and associated signals at other wave- 
lengths is shown in Figure [271 

6. THE ORIGIN OF THE BUBBLE STRUCTURE 

As we have shown in §3l the bilobular Fermi bubble 
structures are apparently well centered on £ = 0°, and 
symmetric with respect to reflection across the Galac- 
tic plane {bottom row in Figure [3]). The bubbles ex- 



tend down to the plane, where they appear even closer 
to £ = 0°. This alignment and north-south symmetry 
are unlikely unless the bubbles originate from the GC, 
and motivate explanations involving Massive Black Hole 
(MBH) activity ( §6.2p or recent starbursts towards the 
GC ( §6.3p . We note that Sco X-1 is approximately cen- 
tered on the north bubble at the present, but this ap- 
pears to be a coincidence. Given the similarity between 
the northern and southern bubbles and the absence of 
any similar feature in the south, we believe the Fermi 
bubbles have a GC origin. 

While the origin of the Fermi bubbles is unknown, a 
rough estimate can be made for their age and the total 
energy required (although the latter quantity depends 
linearly on the gas density in the bubbles, which is poorly 
constrained). From the ROSAT data, we envisage the 
bubbles as hot low density (n ~ 10~^ cm~^) cavities filled 
with ~ 2 keV gas, with (from the Fermi data) height 
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Fig. 21. — The Fermi bubbles at 1 — 5 GeV (the resid ual map o btained by subtracting the SFD dust map and the disk template) 
compared with the WMAP K-band (23 GHz) haze (Do bler fc Finkbeiner 2008 V Top row: The 1-5 GeV map with i = [-45°, 45°] and 
b = [— 90°,0°] (left panel) with the Fermi south bubble edge overplotted in green dashed line (right panel); the same as the left column 
third row panel of Figure [3] Bottom, row: Same sky region as top row but displaying the WMAP haze at 23 GHz (left panel), with the 
Fermi south bubble edge overplotted in green dashed line (right panel). 

pothesis that outflows from a black hole accretion event, 
including possible AGN jets, could form the Fermi bub- 
bles. Section 16.31 focuses on the possibility of a previous 
starburst phase of the Milky Way, with the Fermi bub- 
bles being inflated by the subsequent energetic Galactic 
wind. 



~ 10 kpc, expanding at velocity 'U < 10^ km/s: thus we 
estimate the energy of the Fermi bubbles to be about 
;I^q54-55 g^g^ ^-^^Yi an age of ~ 10'^(v/1000km/s)"^ yr. 

Energetic galactic outflows are common phenomena 
which have b een found in both n earby and high red- 
shift galaxies (jVeilleux et al.l 120051 ). Galactic winds are 
believed to have significant impact on galaxy formation, 
morphology, and their environments. The main sources 
of the energy are stellar winds, super novae, and/or active 
galactic nuclei (AGN); CRs and ma gnetic field pressure 
can also help drive g alactic outflows (jEverett et al.ll2010l : 
Socrates et all l2008h . 

In this section we present some ideas for past activity 
in the GC that could help to generate the shape of the 
Fermi gamma-ray bubbles and the associated signals at 
other wavelengths. Section \6A\ presents and discusses ev- 
idence of past GC activity. Section 16.21 discusses the hy- 



6.1. Observational Evidence of Previous GC Activity 

Observations across the electromagnetic spectrum have 
provided constraints on dynamics and evolution of the 
central gas and stellar populations, and current and prior 
accretion activity of the central MBH. We highlight some 
of the evidence of previous activities towards the GC, 
which may relate to the production of the Fermi bubbles. 

X-ray Reflection Nebulae in the GC — There are indica- 
tions of previous GC activity from X-ray echoes and 
time variability of reflected X-ray lines from cold iron 
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Fig. 22.— Difference maps of the M^MAP haze. Top row: The difference map between the 23 GHz WMAP haze and 33 GHz WMAP 
haze. The ri^'/it panel is the same as the left but overplotted with the south Fermi bubble edge in green dashed line. Bottom row: The 
same as the top row, but showing the difference between the 41 GHz WMAP haze and 33 GHz WMAP haze. We see no apparent structure 
in the difference maps, indicating a consistent spectrum across different (spatial) regions of the haze. 



atoms in molecular clouds around Sgr A* including Sgr 
Bl and B2, Sgr C, and MO 11-0.11 (jSunvaev et al.lll993 : 
Sunvaev fc Churazovlll998l ). The changes in the inten- 
sity, spectrum, and morphology of the fluorescent iron 
nebulae near the GC, observed by ASCA and INTE- 
GRAL are likely due to reflected X-rays from previ- 
ous activity of Sgr A* with high luminosity ^300 yr 
ago. The luminosity is ~1.5xl0^^ erg s~^ in 2 — 200 
keV with a power law spectrum dN/dE (x E~^ with 



7=1.8±0.2 (iPontietal. 



Nobukawa et al.l 



1 )ectru 



Revnivtsev et al.l 12004 : 



200^ 



The changes in the intensities 
and morphologies of hard X-ray nebulosities on parsec 
scales have been discovered (jMuno et al.ll2007l ). 



attributed it to a large-scale bipolar wind, powered by 
a starburst in the GC. In this picture, overpressured 
bubbles rise and expand adiabatically away from the 
injection region in the Galactic plane, driving shocks 
into the surrounding gas by ram pressure. Reflection 
of the shocks increases the density and temperature of 
the post-shocked gas, which has been suggested to have 
T ^ 5 X 10^ K and n ^ 0.1 cm~^ at 2 kpc above the plane 
(jBland-Hawthorn fc Cohenll2003b[ ). Free- free emission in 
the resulting high-temperature plasma produces the ob- 
served X-ray signals. The Galactic Center Lobe (GCL) 
on degree scale has estimated total kinetic energy ~ 10^^ 
erg and a dynamical time scale of ~ 10^ yr. The size, 
energy, and time scales are similar to those of theexpand- 



Outflows: Galactic Cen ter L obe and Expanding Molecu- 
lar Ring - 



ing molecular ring ( EMR) aroun d the GC (JKaifu et al 



S^ ^mnh and lBland-Hawthorn fc Coher] ll272|; | Scoville| | l972 t | Totani| | 2006 |). ISofud (|2Q00bD inter- 



(|2003b[ ) have previously noted the presence of the ex- 
tended bipolar structure in the ROSAT data, and have 



preted the North Polar Spur (NFS) with tens of degrees 
scale to be an outflow from the GC with energy scale of 
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Fig. 23. — The estimated spectrum of IC gamma rays {upper panel) and synchrotron radiation (lower panel) originating from a hard 
electron spectrum along a line of sight 4 kpc above the Galactic center (i.e. b ~ 25°). The steady-state electron spectrum is taken to be a 
power law, dN/dE oc E~^ , with index 7 = 1.5 (solid black), 1.8 (blue dashed), 2.4 (green dotted), and 3.0 (red dash-dotted). In all cases the 
spectrum has a range of [0.1, 1000] GeV. The interstellar radiation field model is taken from GALPROP version 50p, and the magnetic field 
is set to be 10 /xG for the left panel and 5 iiG for the right panel. The data points in the upper panels show the magnitude of the bubble 
emission obtained from template fitting in Figure fT2] ( brown) and Figure [l4] (black) including the "whole bubble" template, as a function 
of energy. The lowest and highest bins contain 3cr upper limits rather than data points with la error bars, due to the large uncertainties 
in the haze amplitude at those energies. For reference, a rectangular cross hatch region shows a approximate spectrum in the same place 
in this and subsequent figures. The data point in the lower panel shows the magnitude of the WMAP haze averaged over b = —20° to 
—30°, for \i\ < 10°, in the 23 GHz K-band (the overall normalizatio n is chosen to fi t this value), and the gray area indicates the range 
of synchrotron spectral indices allowed for the WMAP haze by Dobl er Sz Finkbeinerl (|2008l ). The same population of hard electrons can 
consistently generate the WMAP synchrotron haze and Fermi ICS bubbles. 

- 10^5-56 g^g g^^j ^ij^g g^3g^ig Qf ^ ;lq7 y^^ Ixotanil (|2QQ6h 



suggested tliat all these outflows can be attributed to the 
past high activity of Sgr A* of a duration of ~ 10^ yr, 
comparable to the reasonable estimation of the lifetimes 

of AGNs. 



explained by supernova explosions and the origin of this 
hot plasma, and might be explained as a result of shock 
heating by the wind or AGN activities. Similar diffuse 
hard X-ray emission has b een detected from t he star- 
burst ffala2c^^^l82JStrickl^d^He^^ How- 



Diffuse X-ray Emission — IMuuo et alJ (|2QQ4l ) studied the 
diffuse X-ray emission within ^ 20 pc of the GC in de- 
tail using Chandra observations. The hard component 
plasma with kT ~ 8 keV is spatially uniform and cor- 
related with a softer component with kT ~ 0.8 keV. 
Neither supernova remnants nor W-R/0 stars are ob- 
served to produce thermal plasma hotter than kT ^ 3 
keV. A kT ^ 8 keV plasma would be too hot to be 
confined to the GC and would form a plasma wind. 
A large amount of energy input ~ 10^^ erg/s is re- 
quired to sustain such hot plasma. If the hot plasma 
is truly diffused, the required power is too large to be 



Revnivtsev et al.l (|2009[ ) have resolved ~ 80 percent 



ever, 

of the hard diffuse gas into faint point sources. 

Bipolar Hard X-ray in the GC — Chandra observations 
show the morphology of hot gas with a few keV seems 
to be bipolar with each lobe extending to ~10 pc (see 
Markofj l20ld . and reference therein) . Three explana- 



tions have been suggested for the origin of the lobes: 
thermal wind from the central cluster of massive young 
stars, steady outflows from Sgr A*, o r repeated episodic 
outbursts from Sgr A* (|MarkofJl2010l ) . For the collective 
stellar wind interpretation, it is unclear why the lobes 
only extend up to 10 pc which has a estimated flow time 
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Fig. 24. — The estimated spectrum of IC gamma rays originating 
range, as in the top row of Figure 1231 but with different minimum 
to the data. The three peaks are from ICS of the CMB (left peak) 
ISRF model is taken from GALPROP version 50p. A hard electron CR 
better than a power law electron spectrum extending from 0.1 GeV 

'-«^ 3 X 10^ yr, wliereas tlie star cluster is ~ 6 millon years 
old. Moreover, discrete blobs have been found within 
the lobes, and quasi-continuous winds are hard pressed 
to explain the origin of the blobs. The possibility of 
a transient jet-like feature is intriguing. The jets can be 
produced by accreting the debris of tidal disruption stars 
(see §6.2.11 for more discussion). 

OB Stellar Disk — There are two young star disks that 
have been identified in the central parsec of the GC 
( Paumard et al.ll2QQ6[ ). In situ star formation from dense 
gas accretion disks is favore d over the insp iralings; star 
cluster scenario (see e.g. Bar tko et al.ll2QQ9[ ). The gas 
disks could be formed as a consequence of a large inter- 
stellar cloud captured by the central MBH, which then 
cooled and fragmented to form stars. Interestingly, the 
two star formation events happened near-simultaneously 
about 6 ± 2 My r ago and the two d isks are coeval to 
within ~ 1 Myr ([Paumard et al.ll2QQ6[ ) and little activity 
has occurred since. The two young massive star clusters 
Arches and Quintuplet in the central 50 pc, with similar 
stellar mass, content and mass functions, were formed 



from a hard electron spectrum {dN/dE ex. E~'^) with a limited energy 
and maximum energies. The normalization of the ICS signal is fitted 
, FIR (middle peak), and optical/UV interstellar radiation field. The 
population with a low energy cutoff at about 500 GeV can fit the data 
to 1000 GeV (see Figure [23]). 

( 20061 ) that a global event may cause an increase in the 
rate of star formation, such as a passing satellite galaxy 
that enhanced the clouds' collision frequency and may 
also provide the gas for active star formation. 

The unique characteristics of stellar clusters towards 
the GC have been used to explain the origin of the mag- 
netized nonthermal radio filaments, threads, a nd streaks 
( LaRosa et all 120041 : IVusef-Zadeh et al.l l2004l ) . Collec- 
tive winds of massive W-R and OB stars within such a 
dense stellar environment can produce terminal shocks 
that accelerate relativistic particles. The abundance 
and characteristics of these nonthermal radio filaments 

within the inner 2 degrees of the GC region can be 

I " II 1 

evidence of an earlier s tarbur st (jRosner fc Boddll996f ). 

Yusef-Zadeh fc Konigll (|2004l ) propose a jet model in 



which the characteristics common to both protostellar 
and extragalactic jets are used to explain the origin of 
nonthermal filaments in the Galactic center region. 

6.2. Outflow from Black Hole Accretion Events 
The central MBH in our Milky W ay, with an esti- 



10^ yr ago. It has been suggested in IPaumard et al.l ICillessen et al.l |2009[ ). is currently quiescent, radiating 



mated mass of Mh^/ -4 xlO^ M© fe.g. lGhez et al.ll2008 : 
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Fig. 25. — The estimated spectrum of IC gamma rays originating from a hard electron spectrum, the same as Figure [24l but for the 
500 — 900 GeV energy range. Top row: data points show the separately fitted spectra for the bubble interior and outer shell templates (as 
defined in Figure (TS]), with the template for the disk IC emission given by (left panel) the 0.5 — 1 GeV Fermi map with dust-correlated 
emission subtracted, and (right panel) the simple geometric disk model defined in Figure [S] Bottom row: as top row, but showing the 
separately fitted spectra for the north and south Ferm,i bubbles (as defined in Figure (TS]). 

ing tlie ambient material, producing botli gamma-rays 
and CRs, and appearing more similar to normal low- 
luminosity AGN. If the MBH were radiating at the Ed- 
dington luminosity, it would take only ^ 10^~^ years to 
reach the estimated energy of the Fermi bubbles; for a 
percent level accretion rate (~ 10^^ erg/s), it would take 
~ 10^~^ yr, comparable to the estimated cooling time of 
the electron CRs. 



at about 8 orders of magnitude lower than the Edding- 
ton luminosity {LEdd ^ 10^"^ erg/s). The X-rays of 
Sgr A* is weak and thermal, in contrast to the hard 
nonthermal power law typically observed in most low- 
luminosity active galactic nuclei (LLAGN; see e.g. IHoI 
[2008). Fast X-ray flaring of Sgr A* via nonthermal pro- 
cesses lms_been_dJ^overed by Chandra ACIS observa- 
tion (jBaganoff et al.l |2QQ3[ ). The observed submillime- 
ter/IR bump is seen to flare simultaneously with the 
X-rays. Due to the short timescale of the flares and 
the lack of evide nce of any standard thin accretion disk 
emission (Shakura & Sunvaev 1973) , a magnetic origin 
of the flares has been suggested, either fro m synchrotron 



or synchrotron-self Compton emission (jEck art et al. 
2QQ4I : lYusef-Zadeh et al.ll2QQ6l : iDodds-EdeneTal.. _2QQ9_ ) . 



Both radiative ineflicient accretion flow (RIAF) mod- 

els and outflow-dominat e d models have been invented 

Falcke fc Markofj bOQOl : iBlandford fc BegelmanI Il999 : 



Yuan et al.ll2QQ2[ ). 



Clearly the MBH has not always been so underlumi- 
nous: it may have experienced a long active state in the 
past few million or tens of million years through one or 
more accretion events, driving jets out of the disk, shock- 



6.2.1. Scenarios for MBH Jet Formation in the GC 

Currently, the X-ray flares are the onl y unambiguou s 
AGN-like activity that Sgr A* displays |MarkofJ [2Q1q[ ) . 
However, the synchrotron radio emission is comparable 
to typical LLAGN with compact jets in terms of spectral 
characteristics. It is well known that jets associated with 
accretion disks surrounding black holes can efficiently 
generate high energy particles. BL Lacs can form rel- 
ativistic jets to produce TeV gamma rays, as can mi- 
croquasars. However, we know through multiwavelength 
observations that the central MBH in our Milky Way 
has extraordinarily low bolometric luminosity of ~ 10^^ 
ergs/s, and so is currently in its quiescent dim state, and 
no jets toward the GC have been physically resolved. 
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Fig . 26. — Fermi bubble features in other maps. Top row: The left panel shows the half sky Haslam 408 MHz map (JHaslam et al.l 
Il982l ) with —90° < i < 90°, the middle panel subtracts a simple geometric disk template (shown in the bottom, left panel) to better reveal 
the structures deeper into the Galactic plane. The right panel is the same as the middle panel but overplotted with the Ferm^i bubbles, 
the northern arc, and the Loop I features identified from the 1 — 5 GeV Fermi gamma-ray map (see Figure [4|). The Loop I feature (red 
dotted line) align with the extended diffuse features in the Haslam 408 MHz synchrotron map (known as North Polar Spur). The inner 
and outer edges of the northern arc (dashed blue lines) overlap with two arcs in the Haslam synchrotron map. However, the Fermi bubbles 
have no app arent counterparts in this map. Second row: The same as the top row, but for the Leiden/Argentine/Bonn (LAB) Survey of 
Galactic Hi (JKalberla et al.ll2005l V The middle panel subtracts a simple disk template shown in the bottom middle panel to better reveal 
structures towards the GO. No apparent features have been identified that correlate with the Fermi bubbles and other features in 1 — 5 
GeV Fermi map (there may be some faint filaments morphologically tracing the gamma-ray features). Third row: The same as the top 
row but for the Ha map lF inkbeinei|l2003l ) . The middle panel subtracts a simple disk template shown in the bottom right panel to reveal 
more structures in the inner Galaxy. No corresponding features have been identified morphologically similar to the structures in the 1 — 5 
GeV Fermi gamma-ray maps (color line in the right panel). 
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Detailed examinations of the GCL have shown that the 
gas shell is deep into the disk, and do not support a jet 
origin for that structure (|Lawll2QlQ[ ). However, weak jets 
consistent with the spectrum of Sgr A* ca n be easily hid- 



den b y the blurring of their photosphere (JMarkoff et al. 
2007h . 



The jets related to the energetic accretion in the GC 
can go in any direction; they need not be aligned with the 
minor axis of the Galaxy. This mechanism does not obvi- 
ously provide a natural explanation for the north-south 
symmetry of the Fermi bubbles, and the relatively flat 
gamma-ray intensity inside the bubbles. Also, there is as 
yet no conclusive evidence for the presence of jets toward 
the GC (Muno et al. 2008 ). However, CR rich wide jets 
might have existed in the past, and the Fermi bubbles 
might have been inflated by the jets in a relatively short 
time scale without signiflcant cooling. Although the past 
AGN phase could be the primary heating source of the 
Fermi bubbles, the low density and momentum associ- 
ated with the jets do not readily distribute the thermal 
energy isotropically, making it harder to explain the mor- 
phology and flat intensity of the Fermi bubbles. If a star- 
burst phase in the GC coincided with the energetic jets, 
SNe in the starburst might provide a large injection of 
momentum and turbulence, which could help isotropize 
the energy distribution. 

Accretion of Stars — One way to form jets in the GC is 
for the MBH to tidally disrupt and swallow stars in the 
nuclear star cluster (Hills 19'^ iReeslllQSSJ ). The typi- 
cal picture is that when a star trajectory is sufficiently 
close to the MBH, the star is captured and disrupted 
by tidal forces, and after about an orbital timescale, a 
transient accretion disk forms from the debris of the dis- 
rupted star. The capture rate of the GC MBH has been 
estimated at ~ 4.8 x 10~^yr~^ for main-sequence stars 
and ~ 8.5 x 10~^yr~^ for red giant stars. Sgr A* could 
temporarily behave like an AGN and produce a powerful 
jet by accreting the debris of such stars, which may be 
ejected from surrounding molecular disks. If a 50 solar 
mass star is captured by the MBH in the GC, it gives an 
energy in relativistic protons as high as ^ ]^q54-55 ^^gg 
on a very short timescale (^ 10^~^ yr), at a rate of about 
~ 10^^ ergs/s. 

Accretion of ISM — Quasi-periodic starbursts in the GC 
have been recently suggested as a result of the in- 
teractions betwee n the stellar bar and interstellar gas 
(jStark et al.ll2004l ). In this scenario, gas is driven by bar 
dynamics toward the inner Lindblad resonance, and ac- 
cumulates in a dense gas ring at 150 pc until the critical 
density threshold is reached. Then giant clouds can be 
formed on a short timescale, and move toward the cen- 
ter by dynamical friction. The timescale for this cycle to 
repeat is highly unce rtain but is estimated to be of order 
20 Myr (|Stark et al.l[2004 ). 



Accretion of IMBH — A single 10^ solar mass BH spiral- 
ing in to the GC may als o trigger starbursts and chang e 
the spin of the Sgr A JGualandris fc Merritl (|2009h . 
thereby producing precessing jets. It has been argued 
that one such event happens approximately every 10^ 
years in order to create a core of old stars in the GC, 
of radius 0.1 pc. Chandra has detected more than 2000 
hard X-r ay (2 — 10 keV) p oint sources within 23 pc of 
the GC (jMuno et al.ll2003[ ). Some of them may harbor 
intermediate mass black holes. However, to our knowl- 
edge, no evidence of a collimated outflow has been found 
from the GC. There may be aligned, VLBI-scale (~1 pc 
in the GC) radio knots that have not been discovered. 



Totanil (|2006h 



Models of Past Enhanced Sgr A* Accretion 
suggested a RIAF model of Sgr A* to explain both the 
past high X-ray luminosity and the kinetic luminosity of 
the outflows inferred from observations. The required 
boost factor of the accretion rate ~ 10^~^ with a time 
scale of ~ 10^ yr is naturally expected in the model. The 
induced outflow is energetic enough to support the hot 
{^8 keV) plasma halo towards the GC. The sudden de- 
struction of the accretion flow of Sgr A* is caused by 
the remnant Sgr A East passing through the MBH in 
the pa st ~ 10^ yr. Such a model is claimed by iTotanil 



( 20061 ) to better explain the positron production which 
is required to generate the 511 keV line emission to- 
wards the Galactic bulge obs erved by INTEGRAL/SPI 
( Weidenspointner et al.ll2008[ ). 



Galactic Jets in Starburst Galaxies — In other galaxies, 
powerful AGN radio jets i nteracting strongly with th e 
hot gas have been observed (JMcNamara fc Nulsenll2007l ) . 



In relatively radio-quiet galaxies such as NGC 46 36, 
NGC 708, and NGC 4472, (see e.g. lWang fc m]l2005[ for 
a summary) , faint X-ray "ghost" cavities appear without 
corresponding radio lobes. It has been suggested that 
capture of red giant stars and accretion onto the central 
MBH can power jets/outflows with typical energy ~ 10^^ 
ergs, which is the required energy to form the observed 
cavities. 

Interferometric monitoring of the Seyfert galaxy NGC 
3079 has found evidence of radio jet components un- 
dergoing compression by collision with the clumpy ISM, 
within a few pc of the central engine. This result sup- 
ports the idea that the kpc-scale superbubble originates 
from the spread of the momentum of jets impeded from 
propagating freely. The generalization of this scenario 
provides an explanation for why jets in Seyfert galaxies 
are not able to propa gate to scales of kp c as do jets in 
radio-loud AGN (iMiddelberg et al.ll2007l ). 



Precessing Jets — Several processes may lead to jet 
precession, including magnetic torques, warped discs, 
and gravitational torques in a binary system (see e.g. 
Falceta-Goncalves et al.ll2010[ ). If the momenta of the 
MBH and the accretion disk are not perfectly aligned. 
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the Bardeen-Petters on effect could be a likely mecha- 
nism for precession (jBardeen fc Pettersonlll975[ ). which 
will force the alignment of the disk and the MBH angular 
momentum. 



Irwin fc SeaquistI (Il988f ) suggested that the "bubble" 



structures that have been seen in nearby starburst galax- 
ies could have been blown by a precessing VLBI-scale jet, 
a flat-spectrum radio core was recognized in NGC 3079 
( Baan fc IrwinI Il995l : llrwin fc Sofu3 Il992h . Typically, 
narrow, relativistic and non-precessing jets can carve out 
a hot gas bubble by interaction with ISM and release 
most of their energy far from the GC. Wide jets with 
large opening angles are capable of transferring momen- 
tum into a larger area resulting in the inflation of fat bub- 
bles ( also for precessing AGN jets see [Sternberg fc Soker 
2008h . 



6.2.2. Shocked Shells Driven by AGN Jet 

Relativistic jets dissipate their kinetic energy via in- 
teractions with the ISM. Jets in radio-loud AGN can in- 
flate a bubble composed of decelerated jet matter which 
is often referred to as a cocoon. Initially, the cocoon 
is highly overpressured against the ambient ISM and a 
strong shock is driven into the ambient matter. Then a 
thin shell is formed around the cocoon by the compressed 
medium. The shells are expected to be a promising site 
for particle acceleration. 

As a simple estimate of the dynamics of the expand- 
ing cocoon and shell, we assume the bubbles and shells 
are spherical, and also assume that the ambient mass 
density profile has a form of a power-law given by 
Pa{^) = Pol^/lkpc)""*^"^, where po is the mass density at 
r = 1 kpc. We further assume that the kinetic power of 
the jet, Lj, is constant in time. Under these assumptions, 
the dynamics can be approximately described based on 
the model of stellar wind bubbles. The radius of the 
shock is given by R{t) ~ Qp^ q( LJ2 ij kpc, where 
Po.oi = Po/0-01^ cm~^, L42 = Lj/10^^ ergs s""*^ and 
tj = t/10^ yr. Taking the expected numbers for the 
Fermi bubbles, we get a approximate estimate of the 
bubble size. The total internal energy stored in the shell 
can be expressed as E^ ~ 0.1 Ljt, implying that roughly 
10% of the total energy released by the jet is deposited 
in the shell. 

6.2.3. Buoyant Bubbles 

X-ray images have revealed shock fronts and giant cav- 
ities, some with bipolar structure, in the inner regions 
of clusters, surrounded with X-ray emitting gas. It is 
believed that the power in radio jets is comparable to 
the energy required to suppress c ooling in giant ellipti- 
cal g alaxies or even rich clusters (JMcNamara fc NulsenI 

mm- 

The depressions in the X-ray surface brightness of the 
ROSAT map may themselves indicate the presence of 
empty cavities or bubbles embedded in hot gas, and could 



be interpreted as a signature of previous AGN feedback 
with hot outflows. For adiabatic or supersonic bubbles 
flrst inflated by AGN jets, once the bubble reaches pres- 
sure equilibrium with the surrounding gas, it becomes 
buoyant and rises, because the mass density is lower 
in the bubble than in its surroundings. As the bubble 
moves away from the GC, toward regions with even lower 
density and pressure, it expands. The velocity at which 
the bubbles rise depends on the buoyancy and the drag 
forces. 

The ISM is in turn pushed by the rising bubble, which 
causes a upward displacement behind the bubble called 
"drift". This trailing fluid can give rise to fllaments of 
cool gas. There are indeed filamentary structures in the 
inner Galaxy in both H i and Ha maps. Their identifica- 
tion could support the buoyant bubble scenario. 

6.3. Nuclear Starburst 

Another possible source of dramatic energy injection 
is a powerful starburst in the nucleus. Starburst induced 
Galactic winds are driven by the energy released by su- 
pernova explosions and stellar winds following an intense 
episode of star formation, which create an over-pressured 
cavity of hot gas. The galactic wind fiuid is expected 
to have an initial temperature within the starburst re- 
gion in the range of 10^~^ K even if it has been lightly 
mass loaded with cold ambient gas (jChevalier fc Clegg 
19851 ). The ISM can be swept up by the mechanical en- 
ergy of multiple SN explosions and stellar winds. Large- 
scale galactic outfiows have been observed in starburst 
galaxies both in the local Universe and at high re dshifts 
( Veilleux et aPboOSl : [Bland-Hawthorn et al.lboOTi ). Star- 
burst episodes n ear the GC have been discussed in (e.g. 
Hartmannlll995l ). 



6.3.1. Morphology of Outflows in Starburst Galaxies 

Starburst-driven galactic winds have been studied ex- 
tensively in both multi-wav eband observations and hy- 
drodynamical simu lations (jStrickland fc StevensI I2OOOI : 
Veilleux et aPboOSh . 



AGN or starburst galaxies show bipolar outfiow s 
( Gallimore et al.ll2006l : fsharp fc Bland-Hawthornlboioh . 
The total energy of the superwind has been estimated 
as ^ io^^~^^ ergs, comparable to the estimated energy 
of the Fermi bubbles. The Spitzer Space Telescope has 
found a shell-like, bipolar structure in Centaurus A^ 500 
pc to th e north and south o f the nucleus, in the mid- 
infrared (|Quillen et al.ll2006[ ). The shell has been esti- 
mated to be a few million years old and its mechanical 
energy of 10^^~^^ erg depends on the expanding veloc- 
ity. A small, few-thousand solar mass nuclear burst of 
star formation, or an AGN origin has been proposed to 
explain the formation of the she ll. 



Recently, IWestmoquette et al.l (J2009|) showed that ion- 
ized gas in the starburst core of M82 is dynamically com- 
plex with many overlapping expanding structures located 



37 



at different radii, with compressed, cool, photo-ionized 
gas at the roots of the superwind outflow. Extra-planar 
warm H2 knots and filaments extending more than ~3 
kpc above and below the galac tic plane have also been 
found (e.g. IVeilleux et aIll2QQ9h 

NGC 253 is one of the most famous nearby starburst 
galaxies and is similar to our Milky Way in its overall star 
formation rate, except for a starburst region toward the 
center of the galaxy with spatial extent of a few hun- 
dred pc^^A^_gala£tic_wmd in_J^^^ in 
Ha (JMcCarthvet al.lll987l ). IStrickland fc StevensI (|lOOO) 
discussed the spatial structure in detail in X-ray. The 
wind reaches out to ^9 kpc perpendicular to the disk. 
Filamentary structures as part of projected conical out- 
flow are found in Ha and near-infrared H2 emission. The 
relatively warm gas (~ 10^ K) exists close to the hot gas 
(~ 10^ K). The X-ray filaments tend to be located in in- 
ner regions compared to Ha, and are brighter where the 
Ha emissions are locally weak. The separation between 
the Ha and X-ray filaments is ~70 pc. The spatial distri- 
butions of Ha and X-ray indicate that the inner Galactic 
wind has higher temperature than the outer part. The 
UV e mission seenis to fo rm a shell around the X-ray emis- 
sion (JBauer et al.ll2QQ8[ ). VLBI and VLA observations of 
the nuclear region of NGC 253 at 22 GHz shows no detec- 
tion of any compact continuum source on milliarcsecond 
scales, indicating no low-luminosity AGN in the central 
region of NGC 253. It seems that the starburst region is 
the most plausible explanation f or the source powering 
the wind (jBrunthaler et al.ll2QQ9[ ). 

6.3.2. The Mechanism of Galactic Winds 

The energy and momentum transfer of the Galactic 
wind could be dominated by high thermal and/or ram 
pressure. Materials have been swept-up and entrained as 
part of the wind, and the wind fluid comprises merged 
SN ejecta and massive star stellar wind material with 
ambient gas from the starburst region. Two popular ex- 
trinsic feedback mechanisms have been suggested: ther- 
mally driven winds powered by core-collapse supernovae 
and m omentum- driven winds powered by starburst radi- 
ation (|Coxll2QQ5l ). 

On the other hand, the idea that CRs and magnetic 
fields can help to drive g alactic winds has been kno wn for 
decades (jlpavichl Il975l : iBreitschwerdt et al.l Il99l[ ) . For 
every core-collapse SN, about 10 percent of the energy 
release is converted into CRs {r^ 10^^ erg). These CRs 
interact with the magnetized ISM extensively and ex- 
change momentum through Alfvenic disturbances: the 
characteristic mean free path for a CR proton in the 
starburst phase of the GC is ~1 pc. The effective cross 
section for CR protons and nuclei interacting with the 
ambient gas is much higher than the Thomson cross sec- 
tion for electrons. The luminosity at which CR collisions 
with gas balance gr avity is about 1Q~^ o f the usual Ed- 
dington luminosity (jSocrates et al.ll2QQ8[ ). 



Momentum wind outflowing galactic supershells can 
also be driven b y Lya radiation press ure around star- 
forming regions (jDijkstra fc Loebll2QQ9[ ). The supershell 
velocity can be accelerated to 10^~^ km/s, and it may 
even be able to escape from the host galaxy. The radii are 
predicted to be Tg/i = 0.1 — 10 kpc, with ages tg/^ = 1 — 100 
Myr and energies Egh = 10^^~^^ erg. 

However, the morphology of the galactic wind in 
nearby starburst galaxies inferred from synchrotron. Ha, 
and H i maps is asymmetric about both the galactic mi- 
nor axis and galactic plane, which may suggest the inho- 
mogeneous nature of the ISM. On the other hand, as we 
have shown, the north and south Fermi bubbles are ap- 
proximately symmetric with respect to the galactic plane 
and the minor axis of the disk. The symmetric structure 
of the Fermi bubbles might indicate that they are not 
generated by subsequent interactions with ambient gas 
throughout the wind. 

Furthermore, the typical speed of galactic winds is 
about 200 - 300 km/s. It takes about 5 x 10^ yr for 
CR electrons to reach 10 kpc, but we have not seen 
any evidence of cooling in the gamma ray intensity and 
spectrum. We probably need a faster transport mecha- 
nism of CRs if they were generated from the GC. How- 
ever, the higher the velocity of the wind which entrains 
the CRs, the greater the kinetic energy the wind con- 
tains. The estimated energy of the Fermi bubbles only 
includes thermal energy; if they are actually kinetic en- 
ergy dominated, then the ener gy requirement to form the 
Fermi bubbles is even larger ([Gebauer fc de Boei 12009 : 
Jokipii fc Morfilllll987[ iLerche fc Schlickeiserlll982[ ). 



The estimated supernova r ate in the NGC 253 star - 
burst region is about 0.1 /yr (jEngelbracht et al.lll998[ ). 
Assuming each supernova explosion releases 10^^ erg, and 
10 percent of the energy is transferred to heating the 
Fermi bubbles, this gives a rate of energy injection com- 
parable to ~ 10^^ erg/s. The star formation activity in 
NGC 253 has been underway for about 20 — 30 Myr, and 
is considered to be in a steady state for the CR transport, 



presu mably with a smaller time scale (jEngelbracht et al 
199i). 



It is possible that although the center of our Milky 
Way is currently in its quiet phase, it was recently (in 
the past 10^ yr) in a starburst phase similar to NGC 253 
Heesen et al.l (|2009l ). To our knowledge, however, there 
is no evidence of massive supernova explosions (~ 10"^"^) 
in the past ~ 10^ yr towards the GC, and no apparent 
Galactic wind features have been found in Ha, indicating 
n o strong recent (^ 10^ yr) star formation activity. 

Everett et al.l (|2008l . l2010[ ) compared the synchrotron 



and soft X-ray emission from large-scale galactic wind 
models to ROSAT and Haslam 408 MHz maps. They 
show that a CR and thermally-driven wind could con- 
sistently fit the observations and constrain the launch- 
ing conditions of the wind, including the launching re- 
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gion and magnetic field strength. The comparison of the 
gamma-ray prediction of the wind model with Fermi- 
LAT diffuse emission, especially to the Fermi bubbles 
might have important implications for the CR driven 
wind of our Milky Way. 

6.3.3. Cosmic Rays from a Starburst 

A central starburst might also generate the increased 
population of electron CRs in the GeV-TeV energy range 
required to produce the gamma-ray bubble signals. Star- 
burst galaxies host a greatly increased rate of supernovae 
in the central region. The shocks from supernovae can 
merge and produce energetic galactic scale winds, and 
the enhanced population of supernova remnants is be- 
lieved to accelerate CRs, resulting in orders of magni- 
tude higher CR density than currently expected in the 
GC. It is likely that the GC exhibited comparable CR 
density in the past, with a starburst phase turned on 
by boosted formation of massive stars. As previously 
discussed, during this starburst phase the CR protons 
produced in the inner Galaxy scatter on the ISM with 
very short path lengths, producing gamma rays (via 
TT^ production and decay), and electron (and positron) 
CRs. Although the immediately-produced gamma-rays 
and high density gas and ISRF associated with the star- 
burst phase would not be observable today, the secondary 
electrons might be leftover from the past active star- 
burst phase, and could have been transported to ~10 
kpc by the magnetic field entrained in the Galactic winds. 
GeV and TeV gamma-rays have recently b een detected 
in nearby starburst galaxies NGC 253 (jAcero et al 



2QQ9h . M 82 and Large Magel lanic Cloud (LMC), by 

Fermi-LAT (JAbdo et all l2QlQah. High En ergv Stereo 

scopi c System (H.E.S.S.) (Iltoh et al.lboOTh and VERI 



TAS (JKarlsson fc for the VERITAS collaborationll2QQ9l ) 



support the starburst galaxies as a rich source of high en- 
ergy gamma-rays. If a transient starburst did occur in 
the GC of our Milky Way, and produced a large popula- 
tion of CRs responsible for the observed Fermi bubbles, 
and WMAP haze, what triggered and terminated the 
starburst phase is unclear. 

If the CRs are driven by winds, the halo magnetic 
field can carry CRs along the field lines from the inner 
disk/bulge into the halo, which could help the CR elec- 
trons to reach lOkpc without significant diffusive soften- 
ing of their spectrum. The vertical CR bulk velocity is 
typically hundreds of km/s, which is remarkably constant 
over the entire extent of the disk and for galaxies with 
different mass. In the standard picture, CRs can not 
stream faster than the Alfven speed with respect to the 
static frame of the magnetic field, due to the well-known 



streaming instability (jKulsrud fc Pearcelll969f ). Consid- 
ering nearby starburst galaxies, the vertical CR bulk 
speed has been measured to be vqr = 300±30kms~^ for 
NGC 253. For a typical magnetic field strength of B ^ 
15 /iG and a density of the warm gas of n ^ 0.05 cm~^ 



leads to an Alfven s peed of va = B/^/47Tp ^ 150kms~^ 
(jHeesen et al.ll2009[ ). The super- Alfvenic CR bulk speed 
requires that the CRs and the magnetic field which is 
frozen into the thermal gas of the wind are advectively 
transported together. The measured CR bulk speed is 
the superposition vqr — ^wind + '^a- 

In this picture, the spectrum of electrons injected in 
the starburst region might be harder than elsewhere in 
the Galaxy if the transient starburst phase led to a top- 
heavy mass function of stars. Magnetic turbulence in 
the GC could also be much stronger than in the Solar 
neighborhood: the stronger the turbulence, the faster 
CRs are transported, and the higher the transition en- 
ergy at which synchrotron/ICS losses overtake diffusive 
losses. It has been shown that the high star formation 
rate per area in th e GC leads to short transport times 
(|Becker et al.ll2009l ). 



6.3.4. Constraints on Recent Starburst Activity 

If the Fermi bubbles were generated by previous star- 
burst activity in the GC, we would expect to see many 
more supernova remnants towards the GC than have 
been discovered. Moreover, radioactive ^^Al (half-life 
~7.2xl0^ yr) is believed to be mainly produced by 
massive stars, supernov ae, and novae in the Galaxy 
(|Prantzos fc Diehllll996l ). ^^Al decays into ^^Mg, emit- 
ting a gamma-ray photon of 1808.65 keV. Observations 
of the spectrometer (SPI) on the INTEGRAL gamma- 
ray observatory seem to disfavor the starburst scenario: 
the Galaxy-wide core-collapse supernova rate has been 
estimated at 1.9 d= 1.1 per cen tury from the flux of ^^Al 
gamma-rays (jPiehl et al.ll2006h . Within its half-life, ^^Al 
can only travel ^0.1 kpc with the typical Galactic out- 
flow velocity of ~ 100 km/s. One expects strong ^^Al 
gamma-ray emission concentrated towards the GC, with 
a flux comparable with the total gamma-ray flux from 
the disk, if t he outflow was produced by a starburst 
(JTot anil 120061 ). However, such a strong concentration a t 
the Galactic center is not found (|Prantzos fc Diehllll996[ ). 
indicating that the accretion activity of Sgr A* is more 
plausible as the origin of the mass outflow. Future ob- 
servations of the ratio of the Galactic ^^Fe to ^^Al may 
provide better constraints on the starburst scenario. 

6.4. Other Ways to Generate the Bubbles 

The molecular loops in the GC could possibly be ex- 
plained in terms of a buoyant rise of magnetic loops due 
to the Parker instability. For a differentially rotating, 
magnetically turbulent disk, such magnetic loops can eas- 
ily be formed and rise out of the disk. The typical scale 
of such a loop is Ikpc (JMachida et al.ll2009[ ). 

7. THE ORIGIN OF THE COSMIC RAYS 

It is not necessary that the physical mechanism that 
creates the bubbles also injects the electron CRs respon- 
sible for the Fermi bubbles. It is possible that the bubble 
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structures were formed earlier and the electron CRs were 
injected by an alternative mechanism that then lights up 
the bubble structure with gamma-ray emission. In this 
section, we would like to separate the CR production 
from the bubble formation, and address the spatial ori- 
gin of the CRs. 

In any case, the production mechanism should gener- 
ate electron CRs inside the Fermi bubbles, and also pre- 
vent them from efficiently leaving the bubbles, in order 
to produce the observed "sharp edge". However, elec- 
trons with ~ 10^ GeV diffuse on the order of 1 kpc 
before losing half their energy (jMcQuinn fc Zaldarriagal 



2Q1Q| ). Higher energy electrons lose energy more rapidly 
and so have shorter path lengths; if the gamma-ray emis- 
sion from the bubbles is dominated by IC scattering from 
TeV-scale electrons injected inside the bubble, then the 
sharp edge of the bubbles may be natural. This would in 
turn imply that the gamma rays observed by Fermi are 
largely upscattered CMB photons (starlight and far in- 
frared photons are upscattered to much higher energies), 
which is advantageous for generating such a latitudinally 
extended IC signal. If instead ~ 100 GeV electrons scat- 
tering on starlight are primarily responsible, the electron 
CR density must increase markedly at high Galactic lat- 
itudes to compensate for the falling of the starlight den- 
sity to higher latitude. However, generating a very hard 
CR electron spectrum extending up to (9(TeV) energy 
may be challenging for conventional CR production and 
acceleration mechanisms. These difficulties may be ame- 
liorated if the hard gamma-ray signal from the Fermi 
bubbles is a transient rather than a steady-state solu- 
tion. 

7.1. CRs from the Galactic Center 

As discussed in ^ the electron CRs could be pro- 
duced in the in ner Galaxy by mechanisms such as 
OB associations (jHigdon fc Lingenfelterl |2005|), accre- 



tion events, and supernova explosions, then entrained 
in subsequent jets or outflows and rapidly carried up to 
large scales, avoiding diffusiv e softe ning of the spectrum. 
Breitschwerdt fc Schmutzlerl (|l994l ) suggested that the 
all-sky soft X-ray emission can be explained by delayed 
recombination in a large-scale CR and thermal-gas pres- 
sure driven wind. Such a wind model has been applied to 
the Milky Way which could explain the observed Galactic 
diffuse soft X-ray emission and synchrotron. The model 
indicates that the Milky Way may possess a kpc-scale 
wind. 

The cooling time (denoted r) of TeV-scale electron 
CRs can impose stringent constraints on such models. In 
the Thomson regime, the cooling time scales as (electron 
energy)"^, leading to estimates of r ~ 10^ years for TeV- 
scale electrons. However, where scattering on starlight 
photons is important, electron energies of 100 — 1000 GeV 
are no longer in the Thomson regime, so the scattering 
cross section is suppressed and the cooling time can be 



longer than naively expected. 

Figure [28] shows the variation in cooling time (defined 
as l/{dlnE/dt)) as a function of electron energy and 
height above (or below) the Galactic plane, for the stan- 
dard ISRF model, both with and without the inclusion 
of synchrotron losses (for a simple exponential magnetic 
field profile) ^^. For example, at z = 2 kpc, the cool- 
ing time for TeV electrons is ~ 3 — 4 x 10^ years, rising 
to ~ 7 X 10^ years at z = 5 kpc. Scatterings purely 
on the CMB give an upper bound on the lifetime of 
~ 5 X 10^ (1 TeV/£^) years, but even several kpc from 
the GC, scatterings on the infrared photons dominate 
at TeV-scale energies. The effect of the KN cross-section 
suppression at higher electron energies can be seen in the 
small-z limit where synchrotron losses are neglected, so 
IC scattering of the electrons on starlight is an important 
contribution to the total energy loss rate. 

These relatively short lifetimes, especially close to the 
Galactic plane, may lead to severe difficulties in prop- 
agating CR electrons from the GC out to fill the bub- 
bles. Propagation over such large distances may also 
lead to significant diffusive softening of the electron spec- 
trum, which must be reconciled with the apparent spatial 
uniformity of the bubbles' (gamma-ray) spectral index. 
With electron injection primarily at the GC there is also 
no obvious natural explanation for the flat projected in- 
tensity profile, which seems to require sharp increases in 
the CR density at the bubble walls. 

7.2. CRs from the Bubble Edge 

If the majority of the electron CRs are produced from 
shock acceleration within the edge of the Fermi bubbles, 
the electron CRs in the bubble interior might be left- 
over CRs which undergo cooling after the shock passes 
through. The CRs continue to diffuse inward from the 
shock front while also diffusing outward; if the shock is 
moving faster than the electrons diffuse out, a sharp edge 
in the resulting Fermi bubble gamma-rays is still ex- 
pected. It is also possible that the CRs may be secondary 
electrons, produced by enhanced proton-ISM interaction 
in shocks (within the bubble shell), where protons could 
be ejected from the GC and entrained in the shocks with 
high gas density due to shock compression. 

We can estimate the diffusion path length of 100 — 1000 
GeV electrons, given the lifetimes calculated in §7.11 We 
use the estimated diffusion constant from GALPROP, K = 
5.3 X 10^^ cm^/s at a reference rigidity of 4 GV, and take 
the diffusion coefficient index to be S = 0.43 following 
the results of lAhn et al.l (|2008l ). Then the path length is 
given by. 



E 



ITeV 



0.43/2 



lO^yr 



0.5 



kpc. (4) 



^^ However, note that we have treated the target photon distri- 
bution as isotropic. 
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Fig. 27. — A cartoon picture to summarize the observations of the Fermi bubble structures. Two blue bubbles symmetric to the 
Galactic disk indicate the geometry of the gamma-ray bubbles observed by the Fermi-L AT. Morphologically, we see corresponding features 
in ROSAT soft X-ray maps, shown as green arcs embracing the bubbles. The WMAP haze shares the same edges as the Fermi bubbles 
(the pink egg inside the blue bubbles) with smaller extension in latitude. These related structures may have the same physical origin: past 
AGN activities or a nuclear starburst in the GC (the yellow star). 



Thus we expect the diffusion scale to be smah relative to 
the bubble size, although not negligible. 

Consequently, the electrons in the interior of the Fermi 
bubbles are unlikely to maintain a hard spectrum due to 
diffusion inward from the bubble walls. In this scenario, 
one needs to tune the electron CR distribution to get 
near-flat projected intensity. Although the Fermi bubble 
gamma-rays along any line of sight include contributions 
from both the bubble interior and bubble shell (due to 
the integration along the line of sight), the flat intensity 
and consistent spectrum of the inner bubble and outer 
shell templates (Figure [T6|) implies that the CR spectrum 
in the 3D bubble interior cannot be very different from 
the spectrum at the 3D bubble shell with electron CRs 
generated in situ. 

The edge of the Fermi bubbles might contain MHD 
turbulent fluid with compressed gas and magnetic field. 
Fast magnetic reconnection rather than shocks might 
drive the CR acceleration. Significant magnetic reversals 
within the bubble shell are naturally expec ted, just like 



the heliosheath region o f the solar system (jDrake et al. 



2QQ6I : llnnes et al.l Il997[ ) : the crossing of the termina- 
tion shock by Voyager 1 and 2 may indicate the ac- 
celeration within regions of fast magnetic reconnection 
( Lazarian fc Opherll2QQ9[ ). It is well known that magnetic 



field reversal can cause magnetic reconnection. When 
two magnetic flux tubes of different directions become 
too close, rearrangement of the magnetic field lines takes 
place, converting the energy of the magnetic field into 
energy of plasma motion and heating. Such phenom- 
ena have been investigated extensively in e.g. the solar 
sphere and gamma-ray bursts. In the reconnection re- 
gion, the energetic particles bounce between two mag- 
netic tubes, undergoing first-order F ermi acceleration 
as the dominant acceleration process ([Beck et al.l Il996 : 
Biskamplll986[ lElmegreen fc Scaloll2QQ^ 



On the other hand, maintaining the shape of the bub- 
ble, and preventing it from breaking out, is a no n- trivial 
process. As the bubble rises through the ISM, it tends to 
fatten simply because the fluid moves faster on its sides 
than its front. The classical Ray leigh- Taylor instability 
appears at the leading edge of the bubble, as the inside 
density is lower than that of the ISM; Kelvin-Helmholtz 
instability occurs at the sides of the bubble due to dis- 
continuities in the velocity and density. Although hydro- 
dynamical processes may be capable of stabilizing the 
bubble structures against these instabilities, it has been 
suggested that a large-scale coherent magnetic field could 
help to prevent disruption. We will discuss in §7.31 the 
possibility that the magnetic field can be coherent on the 
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Fig. 28. — The cooling time for electron CRs as a function of energy and height above the Galactic plane, for r 
lines: Synchrotron losses are neglected; equivalent ly, the B-field is assumed to be negligible everywhere. Thick red lines: 
calculation includes synchrotron losses in a magnetic field given by |5| 
GALPROP version 50p, and define the cooling time r = l/{d\nE/dt). 
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: 30e~^/^^P^/iG. We use the standard radiation field model from 



Fermi bubble scale. 

7.3. CR from Diffuse Production in the Bubble 

The hard, spatially uniform spectrum of the Fermi 
bubbles motivates the possibility of some diffuse injec- 
tion of hard CR electrons throughout the bubble volume. 
Such a mechanism, if present, would solve the issues of 
short electron cooling times, relative to the propagation 
time from the GC, and the hardness of the required spec- 
trum. A uniform diffuse injection would give rise to a 
centrally brightened signal, so in models of this type 
there would most likely need to be some other mecha- 
nism increasing the electron intensity at the bubble wall, 
perhaps associated with a shock there. 

This requirement for multiple mechanisms may seem 
unwieldy, and perhaps also unnecessary, if diffusion or 
cooling of electrons produced at the shock can explain the 
flat projected intensity profile and hard spectrum across 
the entire bubble. Given the large size of the bubble, 
however, it is unclear whether production or acceleration 
of electron CRs solely at the bubble walls can give rise 
to a sufficiently centrally bright and hard signal. 

Decay or annihilation of dark matter has previously 
been proposed as a mechanism for diffuse injection of 
very hard electrons and positrons; in particular, dark 
matter annihilation has provided good fits to the data (at 



least in the in ner ^ 15 — 20°) in previous studies of t he 
WMAP haze (iHooper et al.l l2QQ7: 



Harding fc Abazaiian 
2Qld : lLinet al.lboiol : 



2Q1QI: 



Cholis et al.ll2009allb 



Zaldarriaga 



iMcQumn ^.. 

Linden et al.lboiol ). The IC signal 
from electrons produced in dark matter annihilation 
would naively be expected to have approximate radial 
symmetry about the Galactic center and fall off sharply 
at increasing r and z, roughly tracing the distribution 
of dark matter density squared, but smoothed and 
broadened by diffusion of the electrons. This expecta- 
tion is in conflict with the bubble morphology, with a 
gamma-ray distribution elongated perpendicular to the 
Galactic disk, and extending to 10 kpc without much 
change in intensity. However, this naive picture is based 
on isotropic diffusion with a spatially uniform diffusion 
constant, in a tangled magnetic field; including the 
effects of an ordered magnetic field and anisotropic dif- 
fusion can lea d to a much more bubble-like morphology 
(|Doblei<l2010l ). 

Another way to produce CRs in situ inside the bub- 
ble is through magnetic reconnection. Electrons could 
be accelerated directly, or produced as secondaries from 
accelerated protons. The magnetic fields in the under- 
dense bubbles, which may be inflated by AGN outflow, 
may relax into an equilibrium filling the entire volume of 
the bubbles. The timescale depends on the magnetiza- 
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tion and helicity of the outflow and also the properties 
of the ISM. The magnetic field could undergo reconnec- 
tion on a short timescale, converting magnetic energy 
into heat. This mechanism can explain how the bubbles 
could move a large distance through the ISM without 
breaking up. The reconnect ion in the bubbles can also 
accelerate energetic particles, circumventing the problem 
of synchrotron emitters having a shorter lifetime th an the 
age of the bubble they inhabit (|Braithwaitell2QlQ[ ). The 
Fermi bubbles might be initially highly turbulent, with a 
disordered magnetic field far from equilibrium (gas pres- 
sure and Lorentz force are not balanced); the kinetic en- 
ergy would then be dissipated by viscosity, especially in 
the low density bubble. 

The tim e scale of the relax ation to equilibrium can be 
estimated (jBraithwaitdbolOf ) by. 



'^relax ' 



avA 



^4.4 X 10' 







yr, 



(6) 



where va, B and E are the Alfven speed, magnetic field 
and energy at equilibrium. We use the estimated char- 
acteristic values for the Fermi bubbles. The resulting 
time scale is on the order of 10^~^yr. Furthermore, the 
reconnection time scale is orders of magnitude shorter 
than the relaxation time scale: we expect to see ongoing 
reconnection if equilibrium has not already been reached. 

8. POTENTIAL CONNECTIONS TO OTHER OPEN 
QUESTIONS 

8.1. The Cosmic Gamma-ray Background 

Measurement of the in tensity and spectrum of cosmic 



gamma-ray background (|Abdo et al.ll2QlQbl : iFields et al. 



2Qld ) has suggested that instead of rare but intrinsically 
bright active galaxies, it is numerous but individually 
faint normal galaxies that comprise the bulk of the Fermi 
gamma-ray background. This result infers a tighter cor- 
relation between cosmic star formati on history and the 
cosmic gan i ma-ra y background (Thompson et al.l I2QQ7 : 
Lacki et al.ll2QlQl ). Galactic outfiows have been identi- 
fied from near infrared observations at redshift z ^ 2 
( Alexander et al.l |2Q1Q[ ) indicating that such outfiows 
are common features of ultr aluminous infrared galaxies. 
Such outfiows can entrain energetic cosmic rays to escape 
the galaxies, thus providing a way to contribute to the 
cosmic gamma-ray background at energies in the Fermi 
range. The Fermi bubbles may provide local evidence 
for such a scenario. 

8.2. The Origin of Hypervelocity Stars 

Recent surveys of hypervelocity stars (HVSs) have 
found 16 HVSs which are mostly B-type stars with ages 



>(1 - 2) X 10^ yr (see e.g. iBrown et alll2QQ9h . These 
HVSs distributed in the Galactic halo are believed to 
originate fr om the GC involving interactions o f stars with 
the MBH (|Hillslll988l : IYu fc Tremain3l2QQ3l ). Thus en- 
ergetic energy release from the GC which generate the 
Fermi bubbles could b e dynamically re lated to the ejec- 
tion of HVSs. Recent Iv ILu et al.l (|2Q1Q[ ) have shown that 
the spatial distribution of the discovered HVSs is consis- 
tent with being located on two thin disk planes. The ori- 
entation of the planes are consistent with the inner/outer 
clockwise-rotating young stellar (CWS) disk. One pos- 
sibility could be that the HVSs originate from some un- 
known and previously existing disk-like stellar structures. 
HVSs might have been ejected periodically, and related 
accretion events produce jets which generate the Fermi 
bubbles. 

8.3. The Future of the Fermi Bubbles 

What is the future of the Fermi bubbles, are they in 
a "breakout" stage? An interesting possibility is that 
the northern arc and even part of the Loop I feature are 
parts of the rehcs of previous bubbles, and the bubble 
production is a periodic process. The bubbles might be 
fast expanding shocks which might finally expand freely 
into the galactic halo, thus contaminating the ISM with 
entrained hot gas and CRs. An intriguing possibility is 
that CRs and gas released from previous such bubbles 
to the Galactic halo may contribute to the observed dif- 
fuse X-ray and gamma-ray background. In any case, the 
study of Fermi bubbles would have potential implica- 
tions for the understanding of the feedback mechanism 
from the Galaxy. 

8.4. Missing Baryons and High- Velocity Clouds 

N-body/gas dynamical galaxy formation simulations 
have shown that for Milky Way like galaxies, about 70% 
extra baryonic m ass should reside arou nd the galaxy in 
form of hot gas (|Sommer-Larsenl 120061 ). Warm clouds 
are confined by the pressure of the ambient hot halo gas, 
which contains mass at least two orders of magnitude 
more than these warm clouds. The study of the Fermi 
bubbles may provide hints of hot gas feedback from the 
Galaxy, the search for the missing baryons (jBregman 
20071 ) and the puzzle of high- velocity clouds. 



8.5. Metallicity 

Although the Galactic outfiow can not inject a large 
fraction of the ISM, a significant amount of the freshly 
produced metals could be channeled along the galac- 
tic wind. In ordinary photo-ionization it is difficult to 
make the [N II] /Ha ratio exceed about 1.0; shock ioniza- 
tion/excitation is plausible once [N II] /Ha is detected. 
[N II] /Ha has been estimated ~1.5 in NGC 253, less 
than 1.0 in M82 where the ratio trends to increase far 
from the disk. 
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The presence of metals in the IGM has been inter- 
preted as the consequence of energetic metal-rich out- 
flows from galaxies. Active star formation in the in- 
ner Galaxy may contaminates the surrounding ISM, with 
(periodic) Galactic winds entraining the metal-rich gas to 
tens of kpc. The Fermi bubbles may give some hints to 
understanding the feedback and metallicity of the IGM. 
Jets from GC in general do not imply a high metallic- 
ity, and detections of metal rich outflows may essentially 
constrain the energetic injection from jets or Galactic 
outflows from previous starburst towards the GC. 

9. CONCLUSIONS 

We have identified two large gamma-ray bubbles at 
1 < E < 50 GeV in Fermi maps containing 1.6 years of 
data. They have approximately uniform surface bright- 
ness with sharp edges, neither limb brightened nor cen- 
trally brightened, and are nearly symmetric about the 
Galactic plane. The bubbles extend to 50° above and 
below the Galactic center, with a maximum width of 
about 40° in longitude. At |6| < 30°, these ''Fermi bub- 
bles" have a sp atial morphology similar to the WMAP 
microwave haze (|Finkbeinerll2004al : IPobler fc Finkbeiner 
liooa). The RO SAT soft X-ray 1.5 keV map also reveals 
hard-spectrum features that align well with the edges 
of the Fermi bubbles. The similarities of the morphol- 
ogy and hard spectrum strongly suggest that the WMAP 
haze and the Fermi bubbles share a common origin. 

In contrast, the Fermi bubble features are not aligned 
with Loop I or any other feature in the Haslam 408 MHz 
map; while Loop I and other shell structures appear in 
the gamma rays, their spectra are softer than the bubble 
spectrum. Furthermore, there are no convincing features 
spatially correlated with the bubbles in the LAB Hi or 
Ha maps. 

To better reveal the bubble structures, we use spatial 
templates to regress out known emission mechanisms. To 
remove tt^ and bremsstrahlung gammas we use the SFD 
dust map, assuming that the interstellar dust traces gam- 
mas produced by cosmic ray protons (tt^ decay) and elec- 
trons (bremsstrahlung) colliding with the ISM. To trace 
the gamma rays from IC scattering of disk electrons on 
the ISRF, we consider 3 different templates: a simple 
geometric disk model, the 408 MHz Haslam map, and 
the 0.5 — 1 GeV Fermi map after removal of the (dust- 
correlated) TT^ gammas. We verify that our results are 
insensitive to this choice. As an additional cross check, 
we subtract the Fermi diffuse Galactic model from the 
data, finding that this also reveals the bubble structures. 
Although our template fitting technique is subject to sig- 
nificant uncertainties due to uncertain line of sight gas 
and CR distributions, these uncertainties mainly affect 
the intensity profile at low latitudes. For \b\ > 30° and 
1 < £^ < 50 GeV the morphology and spectrum are com- 
pletely consistent for different template choices. Indeed, 
in the 1 — 5 GeV maps, a significant part of the southern 



bubble is easily visible before any template subtractions. 
The Fermi bubbles have an energy spectrum of 
dN/dE ~ E~'^ ^ significantly harder than other gamma- 
ray components. There is no apparent spatial variation 
in the spectrum between the bubble edge and interior, 
and the north and south bubbles have consistent spatial 
and spectral profiles. Both the morphology and spectrum 
are consistent with the two bubbles having the same ori- 
gin and being the IC counterpart to the elec trons which 
generate the microw ave haze seen in WMAP ([Finkbeiner 
2004al : IPobler fc Fi nkbeiner 2008). The spectrum of the 



CR electrons required to generate the Fermi bubbles is 
harder than expected for electrons accelerated in super- 
nova shocks in the disk, and such disk-produced electrons 
would be even softer after several kpc diffusion; the mor- 
phology of the bubble structure is also quite different to 
that of the lower-energy electrons traced by the Haslam 
408 MHz map. Even setting aside the WMAP haze, 
the Fermi bubbles are unlikely to originate from excess 
TT^ emission, as (by construction) they are spatially dis- 
tinct from the SFD dust map, their spectrum is much 
harder than that of the dust-correlated emission, and the 
RO S AT dditdi suggest that the bubbles are hot and un- 
der dense rather than over dense. The morphology of the 
Fermi bubbles, and the overlap of the Fermi bubbles and 
WMAP haze, also strongly disfavor the hypothesis that 
a significant fraction of the high energy gamma rays ob- 
served by Fermi in the bubble region are photons directly 
produced by dark matter annihilation. 

The Fermi bubble structures were likely created by 
some large episode of energy injection in the GC, such as 
a past accretion event onto the central MBH, or a nuclear 
starburst in the last ^10 Myr. We have discussed some 
general possibilities and considerations in this work, and 
found shortcomings in each scenario; it seems likely that 
either significant modifications to one of these ideas, or 
some combination of different mechanisms, will be nec- 
essary. 

Jets originating from AGN activity can potentially ac- 
celerate CR electrons to high energies, and transport 
them rapidly away from the GC; the cooling time of 
electrons at 100 - 1000 GeV is only 10^"^ years, so if 
the CRs are injected and accelerated only in the GC, a 
very fast bulk transport mechanism is required to convey 
them throughout the bubbles before they lose a signifi- 
cant fraction of their energy. However, filling the bub- 
bles completely, with n = 10~^ cm~^ gas, would require 
a mass injection of ~ 10^ solar masses, so in any case 
it is more reasonable for the bulk of the material in the 
bubbles to be swept up and accelerated as the bubbles 
expand. Energetic shocks associated with jets can have 
high Mach number and thus efficiently accelerate CR 
electrons, producing hard spectra with dN/dE ~ £^~^, 
and the total energy required to heat the bubbles is also 
readily achievable by accretion events onto the central 
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MBH. 

However, the north-south symmetry of the bubbles has 
no obvious explanation in the context of an AGN jet: 
there is no reason for one jet to be oriented perpendicu- 
lar to the Galactic plane. The large width and rounded 
shape of the bubbles are also not typical of jets, which 
are generally much more collimated, although a precess- 
ing jet might help explain the wide opening angle of the 
bubbles. If the central MBH becomes active on a rela- 
tively short timescale, the Fermi bubbles may be created 
by a number of past jets, which combine to give rise to 
the symmetric and uniform Fermi bubbles. 

An alternate source for the large required energy injec- 
tion is a nuclear starburst. The wide opening angle of the 
bubbles is not a problem in this case; the bubble shape is 
similar to that observed in NGC 3079, and the X-ray fea- 
tures observed by RO S AT diie similar to those observed 
in other nearby starburst galaxies. However, no corre- 
sponding Ha signal of the Fermi bubbles is observed, in 
contrast to other known starburst galaxies: this prob- 
lem might potentially be resolved if the Ha-emitting gas 
has cooled in the time since the starburst phase (gas hot 
enough to emit the X-rays observed by ROSAT has a 
considerably longer cooling time). Also, generally gas 
filaments and clumps are observed in the X-rays in star- 
burst galaxies, and it would seem that a relic of a past 
starburst should become more clumpy with gas clouds 
and filaments due to cooling of the gas. However, while 
no such structures are obvious in the i^O^^dT maps, the 
signal-to-noise is insufficient to place strong constraints. 

The absence of any such filamentary structures inside 
the Fermi bubbles, on the other hand, argues against 
a hadronic origin for the bubble gamma ray emission. 
Hadronic jets might accelerate protons to high energies, 
and the interactions of these protons with the ISM could 
then produce hard tt^ gammas and secondary e+e~, 
which would scatter on the ISRF to produce more gamma 
rays. In this scenario, however, the gamma ray emission 
should trace the gas density, which we would not expect 
to be smooth and homogeneous. 

Returning to the starburst scenario, the cosmic ray ac- 
celeration in this case would be due to shocks at the edge 
of the bipolar wind. However, the shocks expected in this 
scenario would be relatively weak and slow-moving, and 
thus may not be capable of generating a sufficiently hard 
electron spectrum to reproduce the signal. For example, 
in first-order Fermi acceleration, a shock Mach number 
of ~ 3.3 is needed to obtain an electron spectral index 
of 2.4, as required for the synchrotron explana tion of the 
WMAP haze (see e.g. Ijones fc Ellisonl (|l99lh ). 

The Fermi bubbles have sharp edges, also suggesting 
the presence of a shock at the bubble walls. If the CRs 
producing the gamma rays have a multi-kpc diffusion 
length (which is not expected to be the case for 1 TeV 
electrons, for example), then the edges can still be sharp 



if the bubble edge is moving outward faster than they can 
diffuse. If we assume the Fermi bubbles are projected 
structures from three dimensional symmetric blobs to- 
wards the GC, the flat intensity profile of the bubbles 
requires the emissivity to rise at the bubble walls, but 
remain no n- negligible in the bubble interior; the lack of 
spatial variation in the spectral index may also constrain 
models where the electrons diffuse long distances from 
an injection point. Magnetic reconnection in the interior 
of the bubbles, or some other mechanism such as dark 
matter annihilation, may help maintain a hard spectrum 
throughout the bubbles by accelerating existing lower- 
energy electrons or injecting electrons in situ. 

Dark matter annihilation or decay, while an effective 
mechanism for injecting hard electron CRs at high lati- 
tudes, cannot produce the features in the ROSAT X.-i:d.y 
maps correlated with the bubbles, and would not be ex- 
pected to result in sharp cutoffs in gamma-ray emission 
at the bubble edges. Dark matter annihilation or decay 
may be contributing to the bubbles, or to gamma-ray 
emission in the inner Galaxy that is not well subtracted 
by either the bubble structure template or the models 
for known diffuse emission mechanisms; however, under- 
standing the Fermi bubbles will be a necessary step be- 
fore extracting any such dark matter signal. 

The eROSITA^^ and Planck^^ experiments will pro- 
vide improved measurements of the X-rays and mi- 
crowaves, respectively, associated with the Fermi bub- 
bles, and so may help discriminate between these scenar- 
ios. eROSITA^ which is expected to launch in 2012, will 
provide the first imaging all-sky survey of mid-energy X- 
rays, studying the 0.2 — 12 keV energy range with '-«^ 100 
eV energy resolution and a PSF of 20" . The Planck satel- 
lite, launched in 2009, will greatly improve the measure- 
ments of the WMAP haze spectrum. In addition, AMS- 
02^^ will launch in 2011, and may significantly advance 
our understanding of CR acceleration and propagation, 
and help to refine our interpretation of the Fermi bub- 
bles. 
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